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PREFACE 


This  workshop  will  follow  the  series  which  began  in  Oxford  (2983) 
and  continued  in  Braunschweig  (1984)  and  Trieste  (1987  and  1989). 
Following  the  tradition  of  the  previous  meetings,  the  Workshop  will  be 
devoted  to  recent  advances  in  computational  condensed  matter  physics, 
based  on  realistic  calculations  of  the  electronic  structure  of  polyatomic 
systems. 


The  list  of  topics  to  be  treated  includes:  High  temperature 
superconductivity,  new  materials,  alloy  phase  diagrams,  recent 
developments  in  density  functional  theory,  new  developments  in 
pseudopotential  theory,  approximate  methods.  Car-Parrinello  appraoch  to 
molecular  dynamics,  quantum  Monte  Carlo  treatment  of  correlation  effects, 
time-dependent  correlations  and  excitation  properties,  quantum  dynamics. 
Emphasis  will  be  put  on  both  theoretical  techniques  and  applications. 


As  in  previous  events  in  this  series,  the  Workshop  will  consist  of 
invited  talks,  with  ample  space  devoted  to  discussion.  In  addition,  there 
will  be  one  (or  more)  poster  sesslon(s). 


This  activity  is  co-sponsored  by  the  Scuola  Intemazionale  di  Studi 
Avanzati  (S1SSA,  Tries»p). 


The  Organizers: 

Roberto  Car 
RafTaele  Resta 

Michael  Sehlutcr 


International  Centre  for  Theoretical  Physics  -  Trieste 
Scuola  Internazionale  di  Studi  Avanzati  (SISSA)  -  Trieste 
Consorzio  Interuniversitario  Nazionale  di  fisiea  della  Materia  (INFM) 
Consiglio  Nazionale  delle  Rieerehe  (CNR) 

IBM  -  Italy 

USARDSG  -  European  Research  Office  (London) 


InlmmlioiiaL  Programme  Cflm.mi.tifg,: 

O.K.  Andersen  (Max-PIanek-Institut,  Stuttgart) 

W.  Andreoni  (IBM,  Zurich) 

A.  Baldereschi  (IRRMA,  Lausanne,  and  University  of  Trieste) 

R.  Car  (SISSA,  Trieste) 

S.  Froyen  (Solar  Energy  Research  Institue,  Golden) 

K.  Kune  (Universite  P.  et  M.  Curie,  Paris) 

R.M.  Martin  (University  of  Illinois,  Urbana) 

R.  Needs  (University  of  Cambridge,  UK) 

O.H.  Nielsen  (Nordita,  Copenhagen) 

R.  Resta  (SISSA,  Trieste) 

M.  Scheffler  (Fritz  Haber  Institut,  Berlin) 

M.  Schliiter  (AT&T  Bell  Laboratories,  Murray  Hill) 

E.  Wimmer  (Cray  Research  Institute,  Mendota  Heights) 


International  atomic  Energy  agency 
United  Nations  Educational,  Scientific  and  Cultural  Organiz  ation 

International  Centre  for  Theoretical  Physics 

I.C.T.P.,  P.O,  BOX  586,  34100  TRIESTE,  ITALY,  Cable  CENTRATOM  TRIESTE 


FIFTH  INTERNATIONAL  WORKSHOP  ON 
COMPUTATIONAL  CONDENSED  MATTER  PHYSICS 
(16  -  IS  JANUARY  1991) 


13.00  -  14.30 


L  f  X  C  11  BREAK 
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NEW  SCHEMES  FOR  ELECTRON-ION 

INTERACTION  (Chairperson:  A.  BALDERESCHI) 


14.30  -  14.50 

C-.B,  BA  CHE!  .FT  (llniv.  Trento! :  "Pseudopotentials, 

pseudohamiltomans  and  physical  ions" 

14.50  -  15.00 

Discussion 

15.00  -  15.20 

E.L.  SHTRLEY  (Univ  JUinois.  Urbana):  "GW  quasiparticle 

pseudopotentials 

15.20  -  15.30 

Discussion 

15.30  -  15.45 

BREAK 

15.45  -  16.05 

K  I.AASONFN  (llniv.  Helsinkil:  " Implementation  of  ultrasoft 

pseudopotentials  in  ab-initio  molecular  dynamics" 

16.05  -  16.15 

Discussion 

16.15  -  16.35 

J.M.  SOLER  (Univ.  Autonoma.  Madrid)  "Car-ParrineUo-like 

implementation  of  the  APW  method" 

16.35  -  16.45 

Discussion 

16.45  -  17.00 

COFFEE 

17.00  -  19.00 

POSTER  SESSION 

19.00  -  20.30 

D l N K E R  BREAK 

20.30  -  22.00 

POSTER  SESSION 

THURSDAY,  17  JANUARY 


HIGH  Tc  AND  HIGHLY  CORRELATED 

ELECTRON  SYSTEMS  (Chairperson:  O.K.  ANDERSEN) 


8.45  -  9.15 

E.  TOS  ATTI  (SISSA.  Trieste!:  "Stronelv  correlated  models  for  hieh  T,- 

superconductivity" 

915  -  9.30 

Discussion 

9.30  -  9.50 

W.E.  PICKETT  (NRL.  Washington V  "Lattice  dvnamics  and  structural 

instabihues  in  high  Tc  superconductor^' 

9.50  -  10.00 

Discussion 

10.00  -  10.20 

A. I.  LIECHTENSTEIN  (MP1.  Stuttearii:  "Ab-initio  studs  of  zone- 
boundary  phonons  and  their  coupling  to  electrons  in  YB.0C113O7 
and  Bj  i-x  KxB;03" 

10.20  -  10  30 

Discussion 
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10.30  -  10.50 

S  MASSIDDA  (TRRMA.  Lausaiuiel:  "Electronic  structure  and  positron 

annihilation  m  high  Tc  superconductors:  the  problem  of  the  Fermi 

surface" 

10.50  -  11.00 

Discussion 

11.00  -  11.15 

COFFEE 

11.15  -  11.45 

A  SVANF.  (Ilniv.  Aarhusl:  "Localization  in  the  self- interaction 

corrected  density  functional  formalism' 

11.45  -  12.00 

Discussion 

12.00  -  12.20 

C  A.  B  ALSFTRO  (Bariloche) :  "Static  and  dvnamical  prooerues  of 
quantum  many-body  systems:  Finite  cluster  calculations" 

12.20  -  12.30 

Discussion 

12.30  -  14.00 

L  V K C H  BREAK 

SPECTROSCOPE  OF  HIGHLY  CORRELATED 

FI  FPTRON  SYSTEMS  (Chairperson:  K  KUNC) 

14.00  -  14.20 

H.  ESKES  fl'niv.  Groningen):  "High  energy  spectroscopies  compared 
with  exact  solutions  of  multiband  model  Hamiltonians' 

14.20  -  14.30 

Discussion 

14.30  -  14.50 

M  S  HYBERTSF.N  (AT&T.  Murray  Hill):  "Model  for  low  energy 

electronic  states  probed  by  X-ray  absorption  in  high 
Tc  cuprates" 

14.50  -  15.00 

Discussion 

QUASIPARTICLE  ENERGIES  (Chairperson:  L.  FRITSCHE) 

15.00  -  15.30 

F.  GYGI  (IBM.  Zurich!:  "Self-consistent  electronic  structure  of  a 
vortex  line  in  a  type-II  superconductor' 

15.30  -  15.45 

Discussion 

15.45  -  16.05 

F.  ARYASETIAWAN  (Univ.  Lund):  "  Dielectric  response  and 

quasiparticle  energies  in  Ni” 

16.05  -  16.15 

Discussion 

16.15  -  16.30 

COFFEE 
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NEW  APPLICATIONS  AND  DEVELOPMENTS  OF 
DENSITY  FUNCTIONAL  THEORY  (Chairperson:  P.  ZIESCHE) 

16.30  -  16.50  M.  POSTERNAK  (TRRMA.  Lausanne):  "Quantum-mechanical 

calculations  of  the  spontaneous  polarisation:  Pyroelectric  EiO' 

16.50  -  17.00  Discussion 


17.00  -  17.20  M.  FOULKES  (Imperial  College.  London):  "Alternative  total  energy 

functionals  in  DFT" 


17.20  -  17.30 


Discussion 


17.30  -  17.50 
17.50  -  18.00 


N.  CHK'ITY  (TU  Denmark.  Lvnebvl:  ’’First  principles  energy  density 
and  its  application  to  GaAs  polar  surfaces" 


Discussion 


FRIDAY.  18  JANUARY 

SEMICONDUCTOR  AND  METAL  ALLOYS  (Chairperson:  S.  BARONI1 


8.45  -  9.15 

9.15  -  9.30 


A.  Zl.TNGER  (SERI.  Goldenl:  "First-principle  statistical  mechanics  of 
semiconductor  and  metal  alloys" 


Discussion 


9.30  -  9  50  S.  DE  GIRONCOLT  (IRRMA.  Lausanne):  "Structure  and 

thermodynamics  of  Si-Ge  alloys  from  ab-initio  Mor.te 
Carlo  simulations" 


9.50  -  10.00 
10.00  -  10.15 

10.15  -  10.30 

10.30  -  10.45 

10.45  -  11.00 


Discussion 

K.W.  JACOBSEN  (TIJ  Denmark.  Lvnphvt:  "Alloy  formation  in  the 
effective  medium  theory" 

B.  CHAKRABORTY  (Brandeis  linh^Waliham):  "Ordering  in  Cu-Au 
alloys:  A  new  approach  to  alloy  phase  stability" 


Discussion 

COFFEE 


SURFACES  AND  INTERFACES  (Chairperson:  W,  ANDREON1) 

11.00-  11.30  A.  SF-l.I.ON!  (SISSA.  Triestel:  ''Ab-initio  studies  of  the  Sif  111) 

surface" 

1 1.30  11-5  Discussion 

11.45-  12.05  MD  STIFFS  (NIST  C>a»thrrs\jrc».  "Theory  of  electron  transnussio 

through  interfaces' 


12  05  -  12.15 


Discussion 
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12.15  -  12.25 

O.K.  ANDERSEN  (MPI.  StiittEari):  'Comment  aboui  densitv 

functional  calculations  of  Schottkv  barrier  heights" 

12  25  -  12.45 

R.W.  GODBV  (  Cavendish  Lab,,  Cambridge):  "Exchange  and 

correlauon  m  Schottkv  banners  and  heterojuncnons" 

12.45  •  13.00 

Discussion 

13.00  -  14.30 

LUNCH  BREAK 

14.30  -  15.00 

M  SCHFFFLER  (FHI.  Berlin):  "A  self-consistent  surface  Green- 
f.inction  method,  chemical  trends,  bonding  mechanisms  and 
charge  transfer  of  isolated  adatoms  on  A1  (100)  and  A1  (111)" 

15.00  -  15.15 

Discussion 

DEFECTS  AND  DIFFL'SION  (Chairperson:  J.  BERNHOLC) 

15.15  -  15.35 

P  BLOFCHL  (TBM.  Zurich):  " Ab-initio  calculations  of  self  diffusion 
constants  in  Si' 

15.35  -  15.45 

Discussion 

15.45  -  16.05 

E  SMARGIASSI  tSISSA  Trieste1  'Ab-initio  calculation  of  free 
energies  in  real  materials:  Apphcanon  to  the  vacancy 
formation  in  Si" 

16.05  -  16.15 

Discussion 

16.15  •  16.30 

COFFEE 

LIQUID  AND  , 

AMORPHOUS  SYSTEMS  (Chairrersor:  R.  NEEDS) 

16.30  -  17.00 

D.  HOHL  (NCSA.  Urbans):  "First  pnn  ioles  molecular  dynamics 

simulations  of  amorphous  phosphorus  and  liquid  and  amorphou1 
selenium" 

17.00  -  17.15 

Discussion 

17.15  -  17.35 

G.  SErFFRT  (Tl :.  Dresden):  "Ab-initio  molecular  dynamics 

simulation  of  liquid  Na-Sn  alloy" 

17.35  -  17.45 

Discussion 

CLOSING  SESSION  (Chairperson:  R  NEEDS' 

17.45  -  18.00  V.  HEINE  (Cavendish  Lab.,  Cambridge):  Summarj  and  Outlook 
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ABSTRACTS  OF 


TALKS 


1 1 


Correlated  Ground  States  for  Real  Materials 


D  R  Harr.c-.h 
AT&T  Bell  Laboratories 
Murray  H;"..  New  Jersey  OTST-t,  USA 


With  the  technical  development  of  eiei 
based  on  the  local  density  approximation  at 
appropriate  to  be  exploring  new  paths  tc  greater 
properties,  several  strategies  bsced  on  Monte 
investigated.  These  in:. -is  the  var.ational  appro: 
function  approach  and  the  auxiliary  field  appro 
earn  of  t.tesi  methods,  and  sumtnar!:*  :;s  an, 


ctronic  structure  caicuiar 
today's  advanced  stage 
■  accuracy.  For  ground  s 
Carle  methods  are 
nrh.  the  diffusion  or  Gre 

2  C  .  1  ^  m  *  0  ^  ‘  i  t  j*  w 

t.r/.2ge5.  d :?  a  i  vantages. 


!'M'  e!  aeveioentc!  t .  Sceat.t  pregrto;  ir.  ui.dorrtauc'.r.o  the  aux". 
field  me-:.;;  r.as  increase;  the  probability  ths  :11s  least-envelope;  : e ~ h r. 
may  necome  uselu.  lor  rea:  mat  >r:a‘.«.  Auxiliary  field  prC-.w.' j.,  of  g.. 
state  car;  re  rigorous/-'  fern. ^.ate-i  ££  the  ? n. u * :.«i t.  c-r  m  h'.ff-?!cr.  pc'.a'.'.^v.  A 
ar.alys:*  of  the  ferrr.icr.  s’:g:.  prjc.er:’..  v. h.rh  A  r'.s;cr  barrier  tc 
application  of  auxiliary  f t.-l  s:r.uV;ers.  rtr.  t?  derived  f*orr.  this  formula: 


Based  or.  this  analvsis.  we 


which  ores  r.c  exr.i :  t  the  Cgr  :  •  \  'c  s:._.  gives  »x*e'.>r.t  rtsttlts  for 
tr.erg'.es  and  corrt.ati-.n  t.nc'l.ne  cfsrr.rh  rr  L  .-  ,U-  .mf.-.er  s*»p  trM  ar  > 
in  ate .  f  a  1  ? ,  we  ht-r  r  \Kf 0  e  o  these  sT.m -.atlor.*  •-  hLsca.c-.-t>.?  models 
long  rang*  itre-ictk-.*.  ResuLs.  an.  directions  for  f  u -t  r ?-  development  «1. 
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Calculations  of  the  Properties  of  Solids 
with  Green’s  Function  Monte  Carlo 

David  Ceperley 

University  of  Illinois  at  Urbana-Champaign,  USA 
We  h  ave  recently  calculated  the  equation  of  state  of  solid  Si  with  Green's  Function 
Monte  Carlo.  Pseudo-hamiltonians  were  used  to  eliminate  the  core  degree.1  of 
freedom.  The  cohesive  energy  is  found  to  be  4.66(3)  ev  in  excellent  agreement 
with  experiment  4.62(8)  ev.  This  demonstrates  the  feasibility  of  this  approach  in 
calculating  realistic  systems  and  the  order  of  accuracy  obtainable. 

The  pseudo-Hamiltonian  approach  runs  into  difficulty  with  systems  containing 
transition  metals.  Recently  we  have  performed  Green's  Function  calculation  of 
Si  ,  Sc  and  Cu  atoms  with  a  non-local  pseudopotential  with  the  non-locality 
projected  out  cn  a  highly  optimized  atomic  wavefunction  achieving  reasonable 
results.  However  the  assumed  pseudopotentia!  needs  improvement  to  obtain 
accurate  results. 
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Quantum  Monte  Carlo  Studies  of  Materials: 
Magnetic  Structures  of  BCC  Hydrogen 

Steven  G.  Louie 

Department  of  Physics,  University  of  California 
Materials  Sciences  Division,  Lawrence  Berkeley  Laboratory 
Berkeley,  California  94720 


Abstract 

The  variational  quantum  Monte  Carlo  approach  to  materials  properties 
is  discussed  using  a  recent  application  to  solid  hydrogen  as  example.  Solid 
hydrogen  in  the  atomic  bcc  structure  is  studied  as  a  prototype  single-band 
correlated  electron  system  with  long-range  Coulomb  interactions.  Several 
magnetic  phases  are  investigated  at  densities  between  Wigner-Seitz  radius  of 
rg  *=  1.0  and  rg  =  3.0  a.u.  Correlated  wavofunction  of  the  Jastrow-Slater  form 
with  both  extended  and  localized  single-paritcle  orbitals  are  employed.  A 
transition  from  paramagnetic  to  antiferromagnetic  phase  at  lower  density  is 
found  at  near  rg  *  2.4.  However,  over  the  density  range  studied,  there  is  no 
evidence  for  a  meta6table  partial  ferromagnetic  structure.  If  time  permits, 
results  from  applications  to  other  systems  such  as  D-  centers  in  semiconduc¬ 
tor  quantum  wells  and  in  strong  magnetic  fields  will  also  be  presented. 


dynamical  properties  from  quantum  monte  carlo 

AND  MAXIMUM  ENTROPY 

R.  N.  Silver3,  J.  Gubernati*3,  M.  Jarrell1, 1).  Sins' 

1  Theoretical  Division  and  Manuel  Lujan  Jt  Neuron  Setaertng  Centet  Lei  Alamos  National 
Laboratory.  Lot  Alamos,  New  Mexico  S7545,  USA 
1  Dept,  of  Physics.  University  of  Cinema tti.  Qncmatp.  Ohio  45221,  USA 


We  present  a  new  method1  for  calculating  the  dynamical  properries  of  condensed  matter  sys- 
teres  by  computer  simulations.  It  has  been  successfully  applied  to  the  Anderson  model  for  dilute 
magnetic  alloys.  V.’c  have  obtained  results  for  spectral  functions1,  transport  eoeffecients-  and  dy¬ 
namical  magnetic  susceptibility  which  obey  Kondo  universality  cod  agree  with  experiment  ;  ne 
method  combines  a  quantum  Monts  Curio  method  to  generate  Mautmara  Cneens’  functions  in  imag¬ 
inary  times,  maximum  entropy  and  Hayesian  methods  for  data  analysis  to  analytically  continue  to 
real  rimes,  and  perturbation  theory  to  provide  informative  default  models  needed  as  a  starting  point 
for  the  data  analysis  procedure.  The  method  is  generally  applicable  to  obtaining  dynamics!  tnforma- 
dor.  from  a  wide  rsnge  of  quantum  simulations. 

1  R.  N.  Silver,  ti  al„  Phys.  R **  B41,  23SO  (1990). 

1  R.  N.  Silver,  tt  at.,  Phvj.  Rev.  Lett  65,  496  (1990;. 

*  y.  Jirtll,  rt  Si..  Phys.  Rev  R  Brief  P-tpotta,  January  1991 . 
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Precise  and  Efficient  Ab-initio  Molecular  Dynamics. 
T.A.  ARIAS  M.C.  PAYNE  Cniv.  Cambridge,  UK  and 

J.D.  Joannopoulos  M.I.T.  -  The  Born-Oppenheimer  approx¬ 
imation  divides  the  problem  of  quantum  molecular  dynamics 
into  two  familiar  problems:  (1)  solution  for  the  electronic  wave 
functions  for  a  given  instantaneous  arrangement  of  ions  and 
(2)  the  motion  of  the  atomic  cores  under  the  influence  of  those 
wave  functions.  Combination  of  conjugate-gradient  methods 
to  solve  (1),  standard  molecular  dynamics  to  solve  (2)  and 
an  efficient  method  for  producing  trial  wavefunctions  results 
in  a  scheme  that  is  capable,  with  the  same  computational  ef¬ 
fort  of  traditional  schemes,  of  solving  (1)  three  to  four  orders 
of  magnitude  more  accurately.  The  resulting  ionic  dynamics 
not  only  conserves  energy  over  picosecond  time  scales  with¬ 
out  the  need  for  external  thermostat  controls,  it  also  exhibits 
well  defined  phonon  spectra,  as  determined  from  the  Maxi¬ 
mum  Entropy  Method  (MEM),  in  excellent  agreement  with 
experiment.  Finally,  this  method  may  be  applied  to  systems 
with  large  length  scales  without  complications  involving  con¬ 
vergence  rates  or  stability. 
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Pseudopotentials,  pseudohamiltonians,  and  physical  ions 

G.  B.  Bachelet,  A.  Bosin.  P.  Focher  and  A.  Lastri 

Dipartimento  di  Fisica  and  Laboratorio  di  Fisica  ComputazionaJc 
Universita  degli  Studi  di  Trento,  1-38050  Povo  Trento,  Italy 


The  increasing  interest  in  the  use  of  effective  core-valence  hamiltonians  outside  the 
traditional  context  of  single-particle  schemes  (like  the  Local  Density  Approximation  or  the 
Hartree-Fock  theory)  has  raised  a  number  of  new  questions  in  pseudopotential  theory.  A 
general  question  is  to  which  extent  available  pseudoions,  constructed  within  single-panicle 
schemes,  can  reliably  represent  the  physical  electron-ion  interaction  in  the  context  of  e.g. 
variational  Monte  Carlo,  quantum  Monte  Carlo  or  auxiliary-field  quantum  simulations  of 
realistic  valence-only  systems.  A  second  question,  of  special  interest  for  Quantum  Monte 
Carlo  studies  of  condensed-matter  systems,  is  how  far  one  can  go  in  the  attempt  of  recasting 
contemporary  nonlocal  ionic  pseudopotentials  into  equivalent,  local  pseudo-hamiltonians. 
Theoretical  developments  and  numerical  results  for  the  Periodic  Table  will  be  reviewed. 
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CIV  Quasioarticle  Pseudopotentials.  ERIC  L.  SHIR¬ 
LEY,1  L.  MITAS,  RICHARD  M.  MARTIN,  Univ.  of  Illinois.— 
We  u»e  #gen«tei*ed  version  of  Hedin’s  GW :  to  determine  the 
binding  energy  of  one  electron  outside  an  atomic  core  for  angu¬ 
lar  momenta  of  interest.  We  next  repretent  our  self-energies  by 
semi-local  potentials  from  which  we  generate  pseudopotentials. 
Unlike  traditional  (ingle-body  pseudopotentials,  our  potentials 
systematically  include  core  relaxation  and  core  polarisation  ef¬ 
fects,  as  well  as  “reel’’  core-valence  exchange  and  correlation 
effects.  Tests  of  our  potentials  include  configuration  interac¬ 
tion  and  diffusion  Monte  Carlo  valence  calculations  on  the  Naj 
molecule  and  transition  metal  atoms.  We  also  discuss  using 
“pseudointeractions”  to  improve  pseudopotential  transferabil¬ 
ity. 

1  Supported  by  the  Fannie  and  John  Herts  Foundation. 

5  L.  Hedin  and  S.  Lundquist,  Solid  State  Physict  2S,  1 
(1969),  ed.  F.  Seitz,  D.  Turnbull,  ar.d  H.  Ehrenreich. 


Department  of  Physics 
University  of  Illinois 
1110  W.  Green  Si. 
Urbana,  IL  01601 
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Implementation  of  Ultrasoft  Pseudopotentials 
in  Ab-initio  Molecular  Dynamics 

Kari  Laasonen  and  Roberto  Car 
International  School  for  Advanced  Studies,  S4OI4  Trieste,  Italy 

Changvol  Lee  and  David  Vanderbilt 

Department  of  Physics,  Harvard  University,  Cambridge,  Massachusetts  02138 

A  scheme  for  the  construction  of  ultra-soft  separable  pseudopotentials  recently 
proposed  by  Vanderbilt  is  tested  in  the  context  of  Car-Parrinello  ab-initio  molec¬ 
ular  dynamics  calculations  on  atoms  and  molecules.  For  the  case  of  oxygen,  the 
transferability  of  the  pseudopotential  is  demonstrated  by  comparing  the  calcu- 
lated  properties  of  molecular  0 2  and  O 3  with  those  obtained  from  conventional 
approaches.  Converged  results  are  obtained  using  plane-wave  basis  cutoffs  of  only 
~  25  Ry  dramatically  lower  than  that  required  for  conventional  pseudopoten¬ 
tials.  Forces  can  be  calculated  efficiently,  and  a  molecular-dynamics  simulation 
of  molecular  vibration  is  demonstrated. 

Support  for  this  work  was  provided  by  the  Italian  Cunsiglio  Nazionale  delle 
Ricerche  (CNR)  Grant  No.  89.00006.69,  by  the  European  Research  Office  of  the 
U.S.  Army  Grant  No.  DAJA45-89-C-0025,  and  by  NSF  Grant  No.  DMR-8S- 
172S1.  One  of  us  (D.V.)  also  wishes  to  acknowledge  support  from  the  Alfred  P. 
Sloan  Foundation. 


19 


Car-Parrinello-like  implementation  of  the  Augmented-Piane-Wave  method 
Jose  M.  Soler1  and  Arthur  R.  Williams5 

1Dep.  de  Fisica  de  la  Materia  Condensada,  Univ.  Autonoma.  E-2S049  Madrid.  Spain 

5IBM  Thomas  J.  Watson  Research  Center.  Yorktown  Heights.  New  York  10595 

An  exciting  convergence  of  ideas  is  taking  place  in  the  context  of  effective-single¬ 
particle  electronic-structure  theory.  The  concepts  of  Car  and  Parrinello  are  being  syn¬ 
thesized  with  those  of  Kleiaman  and  Bylander  and  those  of  Slater  to  produce  methods 
capable  of  describing  the  equilibrium  geometries  and  dynamics  of  ’’difficult"  systems, 
such  as  those  with  first-row  and  transition-metal  atoms.  We  have  implemented  one  such 
synthesis  capable  of  very  accurate  description  of  chemically  general  systems.  Molecules 
and  surfaces  can  be  treated  efficiently  without  any  adaptation  and  the  atomic  forces  are 
obtained  directly,  allowing  for  structure  optimization.  Our  work  combines  the  nonlin¬ 
ear  Augmenied-Piaae-Wave  method  with  that  of  Car  and  Parrinello,  but  is  interestingly 
similar  to  that  of  Kleinman  and  Bylander  and  to  recent  proposals  by  Vanderbilt  and 
Blochl.  which  bring  the  pseudopotential  method  to  a  closely  related  form.  These  rela¬ 
tionships  will  be  discussed  and  tve  will  also  present  results  for  several  systems,  small 
molecules,  Cu  phonons,  A1  stacking  faults,  etc. 
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STRONGLY  CORRELATED  ELECTRON  MODELS  FOR  HIGH- 
TEMPERATURE  SUPERCONDUCTORS 

E.  Tosatti 

SISSA/ICTP,  Trieste,  Italy 

I  have  been  asked  to  provide  a  brief  review  of  the  model  Hamiltonians 
which  are  currently  being  considered  in  connection  with  high-Tc  superconductivity. 

Starting  from  the  basic  fact  that  high-Tc  seems  to  arise  by  doping  one 
kind  or  another  of  unconventional  insulator,  I  shall  specialize  to  models  which 
embody  this  feature. 

Issues  of  current  debate,  such  as  absence  or  presence  of  in-plane  hole -hole 
attractive  correlations,  Fermi  liquid  versus  non-Fermi  liquid  behaviour,  and  phase 
separation  will  be  addressed. 
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Lattice  Dynamics  and  Structural  Instabilities  in  High  Tc  Superconductors* 

Warren  E.  Pickett,1  Ronald  E.  Cohen,2  and  Henry  Krakauer3 

1  Cavendish  Laboratory,  Cambridge.  U.K.,  and  Complex  Systems  Theory  Branch,  Naval  Research 

Laboratory,  Washington,  DC  20375-S000,  USA 

2Geophysical  Laboratory,  Washington,  DC  20015-1305,  USA 

3Physios  Department,  College  of  William  and  Mary,  Williamsburg  VA  23185 

A  variety  of  experimental  results  indicate  a  close  relationship  between  the  dynamics  of  the 
lattice  and  the  charge  carriers,  or  of  strong  electron-phonon  coupling,  in  the  high  temperature 
cuprate  superconductors:  Raman  spectra,  phonon  dispersion  curves,  a  low  T  structural 
instability  that  kill  Tc,  phonon  shifts  at  Tc,  changes  in  atomic  structure  at  Tc,  femtosecond 
laser  spectroscopy,  and  unusual  (in  some  cases  large)  oxygen  isotope  shifts,  among  others  In 
an  attempt  to  understand  the  implications  of  electron-phonon  interaction  in  the  cuprates,  we 
have  carried  out  extensive  first  principles,  local  density  functional  calculations  of  phonon 
frequencies,  lattice  instabilities,  and  preliminary  estimates  of  the  strength  and  character  of 
electron-phonon  coupling  in  La2-xMxCu04  and  YBa2Cu3C>7.  The  picture  that  arises  is  one  of 
strong  Madelung-like  shifts  of  potential  due  to  sublattice  displacements  which  lead  to  strong 
ionic  forces  and  can  produce  strong  electron-phonon  coupling.  We  will  also  present 
preliminary  results  on  studies  of  the  low  temperature  tetragonal  (LTT)  structural  transition 
in  Lai  .88®a0.12^'u<-,4  wh'ch  kills  superconductivity  and  produces  a  large  isotope  shift  . 
’Supported  by  the  Office  of  Naval  Research  and  by  National  Science  Foundation 
DM  R-87- 1  9535 . 
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Ab-initio  study  of  zone-boundary  phonons  and 
their  coupling  to  electrons  in  JBagc*sB7  and  ta} _xrxiiOs 

A. Liechtenstein,  C.O. Rodriguez,  O.Jepsen,  O.K. Andersen  and  M.Methfesscl 
Max-PlanckTnstitut  fur  Festkorperforschung,  D-7000  Stuugart  80,  F.R.G. 

The  phonon  frequencies  and  eigenvectors  for  five  S^-modes  in  YBagCugOj  and 
two  R-mode  (breathing  and  tilting)  in  Ba^^K^iOg  (x=0,  0.5)  together  with 
structural  related  properties  were  calculated  from  first-principles  within  the  local 
density  functional  formalism  and  using  the  full -potential  linear  muffin— tin  orbital 
method.  For  YBagCug07  the  full  Sj— dynamical  matrix  was  calculated  using 
63— k  points  in  1/8  Brillouin  zone  for  the  double  cell  and  by  evaluating  the  total 
energy  for  30  displacement  pattern.  We  found  a  large  coupling  between  the 
Op/Og  breathing  mode  in  CuOg- planes  and  the  O ^  mode  in  c-direction.  Normal 
state  phonon  linewidths,  electron-phonon  coupling  constants,  and  softening  of 
phonon  frequencies  due  to  metallic  screening  were  calculated  including 
electron-phonon  matrix  elements.  Furthermore,  comparisson  of  theoretical  and 
experimental  resonant  Raman  scattering  spectra  of  the  5  Ajg  modes  will  be 
presented.  For  Ba ^  _xKiBiOg  at  z=0.5  the  cubic  perovskite  phase  is  stable  while 
for  *=0  an  instability  was  found  for  a  combination  of  breathing  and  tilting  mode 
for  .  The  mechanism  for  this  behaviour  is  explained  in  terms  of  the  competition 
between  the  electrostatic  and  band  energies. 
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Electronic  Structure  and  Positron  Annihilation  in  High  Tc 
Superconductors:  The  Problem  of  the  Fermi  Surface 
S.  Massidda 

Institut  Romand  de  Recherche  Numerique  en  Physique  des  Materiaux  (IRRMAj 
PHB-Ecublens  1015  Lausanne,  (Switzerland) 

The  electronic  structure  of  high  Tc  superconducting  oxides  in  their  normal  state  and  in 
particular  the  existence  of  a  Fermi  surface  (FS)  is  a  crucial  issue  in  these  materials.  While  hov. 
ever  the  agreement  between  band  theory  and  angle-resolved  photemission  spectroscopy  seen- 
to  be  good  for  some  of  these  systems,  still  positron  annihilation  experiments  do  not  show  tin 
agreement  necessary  to  clearly  define  a  Fermi  surface  structure  in  Y Ba^Cn-iOr-i.  We  present 
here  theoretical  calculations  of  the  electron-positron  momentum  density  for  some  high  Tc  su 
perconductors,  based  on  first-principle  LDA  electronic  structure  calculations.  Our  results  for 
Y Ba^Cu^O-!  and  for  the  insulating  parent  compound  Y BaiCn^Ot,  indicate  a  small  overlap  of 
the  positron  wave  function  with  the  CmOj  plane  electrons,  and  consequently  small  signals  du> 
to  the  related  Fermi  surfaces.  By  contrast  we  find,  after  the  folding  of  Vmklapp  terms  a 
cording  to  Lock,  Crisp,  and  West  clear  Fermi  surface  breaks  arising  from  the  Cu  -  0  chair 
bands.  A  comparison  of  the  present  results  with  the  (unfolded)  experimental  two-dimensional 
angular  correlation  of  annihilation  radiation  (2D-ACAR)  recently  measured  in  Geneva  show*,  a: 
overall  agreement  for  )' BajCuzO:  and  a  good  agreement  for  the  insulating  compound,  de-p.v 
the  spurious  LDA  metallic  state.  Our  results  possibly  suggest  the  importance  of  0  vacancie-  ;; 
experiments  performed  on  non-stoichiometrir  Y Ba^Cu-iOi-t  samples.  Predictions  will  be  give; 
for  materials  where  experimental  data  are  still  lacking 


Localization  in  the  Self-Interaction 
Corrected  Density 
Functional  Formalism 


Axel  Svane 

Institute  of  Physics,  Aarhus  University,  Denmark 

and 

Olle  Gunnarsson 

Max-Planck-Institut  FKF,  Stuttgart,  Germany 


A  method  for  performing  fully  self-consistent,  ab-initio  electronic  struc¬ 
ture  calculations  in  the  self-interaction  corrected  (SIC)  local  spin  density 
(LSD)  approximation  is  ,  resen’-d.  The  linear-mufhn-tin  orbitals  method  is 
used  in  the  tight-binding  representation  and  with  the  atomic  spheres  approx¬ 
imation  The  minimization  of  the  SIC  energy  functional  is  achieved  by  the 
steepest  descent  method  with  special  care  taken  to  secure  stabdity  with  re 
spect  to  Unitarian  mixing  of  electron  states  The  method  is  applied  to  the 
3d  monoxides  and  predicts  correctly  that  MnO,  FeO,  CoO,  NiO  and  CuO 
are  antiferromagnetic  insulators  and  that  VO  is  a  nonmagnetic  metal  The 
band  gaps  and  moments  are  significantly  improved  compared  with  the  LSD 
approximation.  The  main  effect  of  the  SIC  is  to  pull  down  the  meta!  d  states 
below  the  oxygen  p-band.  As  a  result,  the  transition  metal  monoxides  emerge 
as  charge-transfer  (rather  than  Mott-Hubbard)  insulators  with  the  highest 
hole-state  having  O  p  character,  however  with  the  symmetry  of  a  metal  d 
6tate.  This  picture  is  in  accordance  with  the  most  recent  interpretation  of 
these  materials  and  at  variance  with  predictions  of  LSD.  Preliminary  results 
for  La2Cu04  are  also  presented. 
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STATIC  AND  DYNAMICAL  PROPERTIES  OF  QUANTUM  MANY  BODY 
SYSTEMS;  FINITE  CLUSTER  CALCULATIONS 

C.A.  Balseiro 
Centro  Atomico  Bariloche 
8400  Bariloche,  Argentina 

Wc  review  a  method  to  calculate  the  ground  state  and  the  frequency 
dependent  correlation  functions  in  finite  many-body  systems  based  on  a  modified 
Lanczos  Method. 

As  examples  of  the  method  we  present  some  results  for  the  electronic 
structure  of  the  Cu02  planes  of  high  Tc  superconductors  and  Green  functions  of 
the  polarized  t-J  model  which  can  be  used  to  study  the  pairing  bags  of  strong 
coupling  superconductors. 
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High  Energy  Spectroscopies  Compared  With  Exact 
Solutions  of  Multiband  Model  Hamiltonians 


Henk  Eskes  and  George  Sawatzky 
Materials  Science  Centre,  Department  of  Solid  State  and 
Applied  Physics,  University  of  Groningen, 
Nijenborgh  18,  9747  AG  Groningen,  The  Netherlands 


ABSTRACT 


In  order  to  extract  information  from  peak  positions  and  intensities  of  high  energy 
spectroscopic  measurements  of  strongly  correlated  mixed  valence  compounds,  it  is  importam 
to  take  explicitly  into  account  the  Coulomb  interactions  as  well  as  the  effects  of 
hybridization.  This  is  the  motive  for  studying  spectra  calculated  by  means  of  an  exact 
solution  (exact  diagonalization  of  small  clusters  or  an  impurity  calculation)  of  multiband 
model  Hamiltonians.  We  will  present  calculations  of  the  photo-  and  inverse  photoemission 
spectra,  as  well  as  O  Is  core  level  spectra  for  the  high  Tc  cuprates.  The  spectra  are  obtained 
for  small  clusters  containing  two  Cu  atoms  (Cu207  and  Cu2Og).  We  discuss  the  sensitivity 
of  the  spectra  to  cluster  size  by  comparing  them  with  results  for  an  (k  dependent)  impurity 
approach  and  results  from  a  Cu04  cluster.  Special  attention  is  paid  to  the  redistribution  of 
spectral  weights  upon  doping. 
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Model  for  Low  Energy  Electronic  States  Probed  by 
X-ray  Absorption  in  High  Tc  Cuprates 

Mark  S.  Hybertsen,  W.M.C.  Foulkes,  and  Michael  Schluter 
AT&T  Bel]  Labors  lories 
Murray  Hill.  NJ.  07974 

E.B.  Siechel 

Condensed  Matter  Theory  Division  1151 
Sandia  National  Laboratories 
Albuquerque,  NM  87185 

X-nty  absorption  near  the  O  K-edge  reveals  important,  chemically  specific,  information  on  the 
electronic  states  near  the  fermi  level.  Recent  bulk  sensative  measurements  by  Chen  and  Sene' 1] 
for  Lai-xSrxCuOt  ahow  the  systematic  evolution  of  the  ipectra  from  x=0,  through  the 
insulator-metal  transition  and  well  into  the  superconducting  regime  (x-0.15).  We  provide  a 
theory  for  these  spectra  starting  from  a  three  band  Hubbard  model  for  the  electronic  states.  The 
parameteri  for  this  model  have  been  previously  derived  from  Local  Density  Functional 
calculations [2],  Direct  diagonalization  studies  show  that  the  model  quantitatively  accounts  for 
the  charge  transfer  gap  and  the  anoferromagnenc  interactions  in  the  parent  insulator[3j.  The 
three-band  Hamiltonian,  including  core-hole  interactions,  is  mapped  into  an  effective  one-bar.d 
model[3].  This  allows  calculation  of  the  X-ray  spectra,  as  a  function  of  carrier  concentration, 
using  exact  diagonalization  techniques  for  clusters  of  adequate  size.  Including  the  role  of  local 
fluctuations  in  carrier  concentration  (e.g.  due  to  fluctuations  in  Sr  concentration)  a  detailed 
comparison  ii  made  to  the  experimental  spectra.  The  resulting  picture  shows  that  the  two  pre- 
edge  peaks  observed  are  precisely  the  0  p-like  states  separated  from  the  Cu  d-like  states  by  the 
charge  transfer  gap  renormalized  by  the  core-hole  excitonic  effect.  The  combined  results  of 
theory  and  experiment  support  the  model  of  the  electronic  states  based  on  a  doped  charge 
transfer  insulator 

[1]  C.T.  Chen,  et  al.,  preprint. 

[2]  M.S.  Hybertsen,  M.  Schluter,  and  N.E.  Christensen,  Phys  Rev.  B  39,  9028  (1989). 

[3]  M.S.  Hybertsen,  E.B.  Stechel,  M.  Schluter,  and  D.R.  Jennison,  Phys.  Rev.  B  41,  11068 

(1990,. 


Self-Consistent  Electronic  Structure  of  a 
Vortex  Line  in  a  Typc-II  Superconductor 

Francois  Gygi*  and  Michael  Schluier 

AT&T  Bell  Laboratories 
600  Mountain  Avenue 
Murray  Hill,  NJ  07974 


The  electronic  structure  of  a  vortex  line  in  a  type-11  superconductor  is  calculated  self- 
consistently  in  the  framework  of  the  Bogoliubov-de  Gennes  theory.  The  Debye 
frequency,  the  Fermi  velocity  and  the  coupling  constant  of  the  electron-electron  attractive 
interaction  are  used  as  microscopic  input  parameters.  The  resulting  quasiparricle 
excitation  spectrum,  the  pair  potential  and  the  current  distribution  are  studied  as  function 
of  temperature,  and  can  be  used  to  define  a  coherence  length  and  to  determine  the 
magnetic  penetration  depth.  The  local  density  of  one-particle  excitations,  calculated 
from  the  quasipanicle  amplitudes  explains  the  results  of  Scanning  Tunneling  Microscopy 
(STM)  experiments  by  Hess  et  al.  JPhys.  Rev.  Lett.  62,  214  (1989)1  on  KbSci  The  main 
spectroscopic  features  in  the  experimental  results  are  caused  by  bound  states  in  the  vortex 
cores.  Spatial  distortions  of  the  bound  state  wavefunctions  caused  by  neighbouring 
vortices  and  by  the  crystallline  lattice  are  discussed  in  terms  of  a  simplified  rwo-bar.d 
model.  In  the  case  of  NbSe2,  the  resulting  local  density  of  states  has  a  characteristic 
"star"  shape  in  real  space,  whose  orientation  is  energy-dependent,  in  agreement  with 
recent  STM  experiments  [Phys.  Rev.  Lett.  64, 2711  (1990)]. 


Present  address:  IBM  Research  Division,  CH-8803  Ruschiikor.  Switzerland 
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Dielectric  response  and  quasiparticle  energies  in  Nickel 
Ferdi  Arvasetiawan  and  Ulf  von  Barth 
Lund  University,  Solvegatan  14A,  S-22362  Lund,  Sweden 

Field-theoretic  methods  in  the  many-body  theory  of  solids  was  an  important  area  of  research  in 
the  early  sixties.  Unfortunately,  neither  the  computational  facilities  nor  the  numerical  techniques 
of  those  days  were  adequate  for  working  out,  in  real  systems,  the  predictions  of  even  the  simpler 
approximations  proposed  at  the  time.  Over  the  past  fifteen  years  the  situation  has  improved 
considerably.  The  dielectric  response  and  the  connected  optical  absorption  spectrum  of  elemental 
semiconductors  has  been  worked  out,  at  least  approximately,  in  the  well  known  RPA  meaning 
free  particles  responding  to  the  applied  plus  the  induced  fields  [1,2,3].  Simple  attempts  to  go 
beyond  the  RPA  and  include  also  the  particle-hole  interaction  have  also  been  made  [3,4].  The 
results  are  quite  encouraging  compared  to  those  obtained  from  the  standard  methods  based  on 
one-particle  band  theory  and  a  Golden-Rule  expression  for  non-interacting  particles.  Recently, 
several  people  [5,6,7]  have  used  the  RPA  response  of  silicon  in  order  to  calculate  the  one-electron 
excitation  spectrum,  i.e.  the  quasiparticle  band  structure,  of  the  system  in  an  approximation 
proposed  by  Hedin  (8]  in  the  mid  sixties  and  often  referred  to  as  the  GW  approximation.  Again 
the  results  are  superior  to  those  of  standard  band  theory  methods.  For  one  thing,  the  calculated 
band  gaps  are  close  to  experiment. 

In  the  present  work  we  have  concentrated  our  efforts  on  transition  metals  using  nickel  as 
the  prototype  case.  As  our  starting  point  for  the  many-body  perturbation  expansion  we  have 
chosen  a  zero:th-order  Green’s  function  obtained  from  a  self-consistent  band-structure  calcula¬ 
tion  based  on  the  local-density  approximation  (LDA).  This  Green's  function  is  known  to  give 
accurate  charge  and  spin  densities  although  the  quality  of  the  one-electron  energies  are  relatively 
poor.  The  d-band  width  is,  e.g.,  almost  30a  factor  two  too  large.  We  proceeded  to  calculate 
the  dielectric  matrix  and  the  closely  related  screened  Coulomb  interaction  within  the  RPA.  Due 
to  consistency  requirements,  as  discussed  by  Kadanoff  and  Bavm  [9,10],  we  have  so  far  made 
no  attempts  to  include  the  effects  of  particle-hole  interaction  either  in  the  response  function  or 
in  the  self-energy.  FVom  the  response  function  (essentially  the  dielectric  matrix)  we  obtain  the 
energy-loss  spectrum,  the  optical  absorption  spectrum  and  the  closely  related  optical  conduc¬ 
tivity  spectrum.  Contrary  to  some  suggestions  we  find  that  the  local-field  effects  are  small  - 
especially  at  small  frequencies  and  small  momentum  transfers.  Therefore,  our  optical  conduc¬ 
tivity  spectrum  closely  reflect  structure  in  the  one-electron  spectrum  of  the  starting  Green’s 
function  and,  consequently,  the  conductivity  peaks  occur  at  different  energies  compared  to  their 
experimental  counterparts  [11].  The  discrepancies  are  of  the  order  of  1.5  eV  .  Starting  from 
a  Green’s  function  constructed  from  LD  wave  functions  but  instead  using  either  experimental 
one-electron  excitation  energies  or  our  calculated  quasi-particle  energies  (see  below)  would  result 
in  an  optical  absorption  spectrum  in  much  closer  agreement  with  experiment.  Such  a  procedure 
would  ,  however,  require  further  theoretical  justification. 

FVom  the  calculated  screened  interaction  and  the  LDA  Green’s  we  computed  the  GW  ap¬ 
proximation  for  the  electronic  eelf -energy.  FVom  this  we  can,  e.g_  obtain  the  spectral  function 
of  the  interacting  one-electron  Green’s  function,  the  quasi  particle  energies  i.e.  the  positions  of 
well  defined  peaks  in  the  spectral  function  as  well  as  extensive  quantities  like  the  charge  density 
and  the  total  energy  per  atom  of  the  solid.  The  resulting  quasi  particle  energies  are  in  general 
closer  to  the  Fermi  level  as  compared  to  the  eigen  energies  of  the  LDA.  The  deeper  d  states 
move  upw  ards  by  0.5  to  1.0  eV  resulting  in  an  approx.  25are  within  a  few  tenths  of  an  eV  from 
the  experimental  results  :  12.13.1 4  . 

The  differences  between  corresponding  quasi  particle  energies  of  the  spit. -up  and  the-  spit, 
down  channels  are  reduced  (i0-3uLDA  but  we  are  still  far  from  the  SOexperiment.  Neither  do 
our  spectral  function  show  a  satellite  peak  corresponding  to  the  famous  two  iioie  state  15. 10 
6  eV  below  the  Fermi  level  as  indicated  by  experiment  or  by  theoretical  calculations  on  model 
systems.  We  believe  we  understand  these  shortcomings  and  our  results  demonstrate  the  necessity 
to  go  beyond  the  GW  approximation  for  certain  properties  in  the  narrow  band  systems. 
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physically  inspired  approximation  schemes.  We  believe  this  knowledge  to  be  instrumental  in 
attacking  more  complicated  but  also  more  rewarding  systems  such  as  e.g.  nickei  oxide. 
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Quantum-mechanical  Calculation  of  the  Spontaneous  Polarization:  Pyroelectric  BeO 
M.  Poitemak,  R.  Rests',  and  A.  Baldereschi 
Institut  Romand  de  Recherche  Numerique  en  Physique  des  Materiaux  (IRRMA), 
PHB-Ecublens,  CH-1015  Lausanne,  Switzerland 

Crystals  which  can  carry  a  spontaneous  macroscopic  polarization  constitute  a  special  class 
of  useful  solid-state  materials.  Their  electrical  polarization  P  can  be  modified  by  changing  the 
temperature,  as  in  pyroelectrics,  or  by  applying  external  electric  fields,  as  in  ferroelectrics  where 
even  the  direction  of  the  polarization  may  be  altered.  Only  variations  of  P  are  experimentally 
detectable,  through  measurements  of  pyroelectricity  or  feiroelectricity.  The  key  idea  of  our  ap¬ 
proach,  in  order  to  calculate  P,  is  to  make  use  of  the  differential  concept  also  in  the  theoretical 
treatment  [lj.  Two  steps  ate  necessary,  which  constitute  in  fact  an  alternative  definition  of  spon¬ 
taneous  polarization:  (i)  the  calculation  of  the  difference  AP  between  two  phases  of  the  came 
material,  and  (ii)  a  suitable  symmetry  argument.  The  simplest  material  where  bulk  polarization 
is  allowed  by  symme.ry  is  BeO:  it  has  the  hexagonal  wurtzite  structure  and  its  constituents  are 
first-row  elements.  The  crystal  has  a  preferred  axis,  which  is  the  direction  of  spontaneous  polar 
ization.  This  latter  quantity  is  obtained,  using  all-electron  local-density  FLAPW  calculations 
and  the  two  steps  outlined  above.  Resulting  values  of  P  are  discussed,  and  compared  to  existing 
experimental  data. 

*  Permanent  address:  Scuola  Intemazionale  Superiore  di  Studi  Avanzati,  Trieste.  Italy. 

[1]  M.  Postemak,  A.  Baldereschi,  A.  Catellani,  and  R.  Resta,  Phvs.  Rev.  Lett.  64,  1777  (1990; 
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ALTERNATIVE  TOTAL  ENERGY  FUNCTIONALS  IN  DENSITY 
FUNCTIONAL  THEORY 

W. Matthew  C.  Foulkes 
Blackett  Laboratory, 

Imperial  College  of  Science,  Technology  and  Medicine 
Prince  Consort  Road,  London  SW7  2BZ,  England 

It  has  recently  been  realized  that  the  Hohenberg-Kohn  energy  functional, 
which  gives  the  total  energy  of  the  interacting  electrons  in  an  atom,  solid  ot 
molecule,  as  a  functional  of  the  electron  density,  is  only  one  of  several  possible 
alternatives.  In  self-consistent  calculations,  it  may  sometimes  be  computationally 
convenient  to  use  one  of  these  alternative  functionals  in  place  of  the  original;  and  in 
non-selfconsistent  calculations  (when  one  is  trying  to  obtain  quick  but  reliable 
results  based  on  some  input  guess  at  the  potential  or  electron  density)  it  may  also  be 
more  accurate.  This  talk  examines  the  freedom  available  in  choosing  alternative 
energy  functionals,  describes  some  of  the  ways  in  which  they  have  been  used  up 
until  now,  and  makes  some  suggestions  as  to  how  they  may  be  used  in  future. 
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First  Principles  Energy  Density  and  its  Applications  to  GaAs 
Polar  Surfaces,  N.  Chetty*  and  Richard  M.  Martin,  University  of  Wi- 
aoii  at  Utbana-Chtmpaign,  Department  of  Physics  and  Materials  Research 
Laboratory,  1110  W. Green  Street,  Urbana,  Illinois  61801 

Density  Functional  Theory  has  been  successfully  applied  in  solid  state  cal¬ 
culations  where  is  has  been  customary  to  extract  a  single  number,  the  total 
energy.  Within  this  formalism,  we  show  that  one  can  in  addition  define  a  local 
energy  density  £(r)  which  is  a  sum  of  kinetic,  Maxwell,  exchange-correlation 
and  other  contributions  that  depend  on  the  external  potential.  Despite  the 
nonuniqueness  of  £(r),  we  show  that  it  can  be  used  to  calculate  well-defined 
energies  from  integrals  over  regions,  e.g.  surfaces.  We  present  explicit  calcu¬ 
lations  for  GaAs  (100)  surfaces  where  we  demonstrate  the  efficiency  of  using 
£(r)  to  determine  the  energies  of  the  two  inequivalent  surfaces  from  a  single 
csuculation.  We  also  calculate  the  formation  enthalpy  of  single  isolated  GaAs 
(111)  and  (ITT)  surface  which  is  not  possible  using  conventional  total  energy 
methods.  We  are  able  to  address,  for  the  first  time,  questions  related  to  the 
itability  of  GaAs  surfaces. 

•  present  address,  Laboratory  of  Applied  Physics,  Technical  University  of  Den¬ 
mark,  Building  307,  DK-2800  Lyngby,  Denmark 
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First-principle  statistical  mechanics  of 
semiconductor  and  metal  alloys 


A.  ZUNGER 

Solar  Energy  Research  Institute 
1617  Cole  Boulevard  16/3,  Golden  CO  80401,  L.S.A. 


Ordering  ia  Cu-Au  Alloy*:  A  New  Approach  to  Alloy  Phage  Stability 
Bulbul  Chakraborty,  Phytic*  Department,  Brandeit  Univertity,  Waltham,  MA  02254-0110 

Abetract 

New  raeult*  concerning  ordering  in  Cu-Au  alloy*  will  be  presented.  The  ordering  hu  been  etudied 
ueing  a  formaliem  baaed  on  the  effective  medium  theory  (EMT)  of  chemical  binding.  The  EMT 
i*  an  approximate  method  baaed  on  density  functional  theory.  It  i*  a  variational  approach  to 
calculating  total  energiee  of  arbitrary  atomic  configuration*  and  becauee  of  ita  limplicity  it  can 
he  need  efficiently  in  numerical  limulttioni.  It  will  be  ihown  that  both  ground  etati  propertiae 
and  order-dieorder  traneition  Umperaturee  of  the  Cu-Au  alloy*  com*  out  in  very  good  agreement 
with  experiment.  Raeulta  of  Monte  Carlo  eimulation*  of  the  phaee  diagram  will  be  presented  and 
the  origin!  of  the  ordering  traneition*  will  be  diseueeed.  Finally  the  virtue*  and  drawback*  of  EMT 
will  be  diecueeed  in  the  context  of  exiating  ab-initio  density  functional  calculation  of  alloy*. 
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Ab-initio  studies  of  the  Si(lll)  surface 


A.  Selloni 

International  School  for  Advanced  Studies,  1-34014  Trieste,  Italy 

A  number  of  applications  of  the  ab-mitio  molecular  dynamics  scheme  to 
semiconductor  surfaces  are  presented.  We  shall  focus  in  particular  on  some 
properties  of  the  defect-free  and  defective  clean  Si(lll)  surface,  including: 

(i)  The  2x1  reconstruction  of  Si(lll),  and  in  particular  the  transformation 
process  from  the  bulk- terminated  lxl  to  the  reconstructed  2x1  surface,  and 
the  details  of  the  resulting  structure. 

(ii)  The  influence  of  surface  point -defects,  namely  surface  vacancies,  on  this 
transformation  process. 

(iii)  Surface  vacancies  and  their  connection  to  the  recently  discovered  v  3  x  %  3 
reconstruction  of  the  clean  Si(lll)  surface. 

(iv)  The  structural  and  electronic  properties  of  an  isolated  surface  vacancy 
on  Si(  1 1 1  )2x  1. 

Our  studies,  besides  demonstrating  the  feasibility  of  ab-initio  calculations  for 
rather  complex  systems,  such  as  surface  defects,  bring  out  several  interesting 
specific  features  which  will  be  discussed  in  detail  in  this  talk. 

This  work  was  done  in  collaboration  with  F.  Andlotto,  W.  Andreoni.  R. 
Car.  M.  Parrinello  and  E.  Tosatti.  It  was  supported  in  part  by  the  Italian 
Ministry  of  the  University  and  Scientific  Research  under  the  SISSA-CINECA 
collaborative  project,  by  the  Italian  CNR  under  grant  no  89.00006.69.  by  the 
National  Institute  for  the  Physics  of  Matter  (INFM),  and  by  the  European 
Research  Office  of  the  US  Army. 
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Theory  of  Electron  Transmission  through  Interfaces.  M.  D.  Stiies,  National  In¬ 
stitute  of  Standards  and  Technology  and  D.  R.  Hamann.  AT&iT  Bell  Laboratories 
As  material  growth  techniques  improve,  the  detailed  scattering  properties  of  atomical.;, 
abrupt  interfaces  become  more  and  more  important.  For  many  of  these  interfaces  simple 
methods  to  calculate  scattering  properties,  like  the  effective  mass  approximation,  car.  not 
give  reliable  results.  We  have  developed1  a  method  of  calculating  electron  transmission 
and  reflection  probabilities  of  such  interfaces  based  on  a  layer-by-layer  construction  of 
the  interfacial  electronic  states.  The  method  is  based  on  a  linearized-augmented-piar.e- 
wave  implementation  of  the  local-density  approximation.  The  electronic  potentials  fcr 
the  layers  are  taken  from  separate  self-consistent  calculations  assuming  that  each  layer 
is  either  bulk-like  or  is  similar  to  a  section  of  a  supercell.  Using  these  potentials,  varia¬ 
tional  electronic  states  are  constructed  in  each  iayer,  and  joined  by  matching  the  values 
and  slopes  at  the  boundaries  between  adjacent  layers.  The  scattering  states,  and  hence 
the  transmission  and  reflection  probabilities,  are  found  by  expanding  the  states  in  trie 
bulk-like  layers  in  terms  of  the  appropriate  generalized  Bloch  states.  The  transmission 
properties  of  NiSi2/Si(lll)  interfaces  are  of  particular  interest  because  not  only  can  good 
interfaces  be  grown,  but  interfaces  with  two  different  relative  orientations  can  be  grown. 
Our  prediction,2  of  substantial  differences  in  transmission  for  the  two  interfaces  has  mo¬ 
tivated  several  groups  to  look  for  them. 

1.  M.  D.  Stiles  and  D.  R.  Hamann,  Phvs.  Rev.  B  38,  2021  ( 1 9SS ) . 

2.  M.  D.  Stiles  and  D.  R.  Hamann,  Phvs.  Rev.  B  40,  1349  (19S9i. 
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Comment  about  density  functional  calculations  of 
Schottky  barrier  heights 


O.K.  ANDERSEN 

Max-Planck-Institut  fur  Festkorperforschung 
Heisenbergstrasse  1,  D-7000  Stuttgart,  Germany 
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EXCHANGE  AND  CORRELATION 
IN  SCHOTTKY  BARRIERS  AND  HETEROJUNCTIONS1 

R.W.  Godby,  Cavendish  Laboratory,  University  of  Cambridge 
Madingley  Road,  Cambndge  CB3  OHE,  U.K.. 

Exchange  and  correlation  at  semiconductor  interfaces  is  discussed,  with  particular 
reference  to  the  contrast  between  the  local  quasiparticle  energies  and  local  density-function.:! 
theory  (and  LD A)  eigenvalues.  This  includes 

•  Consideration  of  the  exchange-correlation  potential  in  the  neighbourhood  of  such  an 
interface.  Slow  variation  of  the  potential  within  the  depletion  layer  can  give  nse  to  bending  of 
the  DFT  bands  relative  to  the  quasipanide  bands2. 

•  A  computational  study  of  many-body  effects  in  a  Al-GaAs(l  10)  Schottky  barrier.  The  G'U 
calculation  shows  significant  effects  on  the  semiconductor's  excitation  energies  in  the  vicinity  of 
the  metal  which  are  absent  from  conventional  local-density  approximation  (LDA)  calculations'. 

1  Work  done  with  J.  Charlesworth  and  RJ.  Needs  (University  of  Cambridge);  LJ.  Sham  (University  of 
California  San  Diego)  and  M.  S'hliiter  (AT&T  Bell  Laboratories). 

2  R.W  Godby,  L.J.  Sham  and  M.  Schliiter,  Ptiys  Rev.  Leu  65.  2083  (1990). 

2  J.P.A.  Charlesworth,  A.  Oschlies,  R.W’.  Godby.  R.J.  Needs  and  L.J.  Sham,  to  appear  in  Proceedings  o'  the 
20ih  International  Conferen '•  on  the  Physics  of  Semiconductors.  World  Scienufic  Publishers.  1990.  J.P.A 
Charlesworth,  R.W.  Godby  and  RJ.  Needs,  to  be  published. 


A  self-consistent  surface  Grcen-functiori  method; 
chemical  trends,  bonding  mechanisms  and  charge  transfer  of  isolated 
adatoms  on  Al  (100)  and  A1  (111) 

M.  SCH HITLER 

FVitz-Haber-Institut  der  MPG,  Berlin,  Germany 

We  describe  the  bar.:  aspects  of  a  new,  leif-vOnsisTent  Greer.-  function  method  v.  _ 
allows  to  calculate  the  dcnsuy  of  states,  election  density,  anti  related  quoitUf.es  to 
localized  perturbation  •«.*;.  an  isolate  !  adsorbate  or  ar  intfusic  surface  bi-fett.  a 
crystal  surface.  The  atomistic  nature  and  the  macroscopic  extension  of  the  cubit: 
arc  fully  taken  into  account.  The  method  is  bail'd  or.  the  density-functional  theory 
combines  seveia.  iceas  from  recent  theoretical  ccve.opmer.ts.  a-'  for  e.sampic  front 
layer-KKR  Green  f  unction  method,  fro:::  ab-:<utiv  pseud  ^potential  theory,  and  .: 
the  self-consistent  defect  Grcen-f  iu.. .or.  me: nod 

in  particular  «:  d-scuss  apnlira' .on-  of  th.e  I’-.eth-xi.  namely  the  acsorpt.oii  of  or 
Ill  adatoms  No  —  C .  ot:  A!  1.1;  a::u  Al  GCo  .  V\r  ade.'.--i  c  rrrcr.t  eon  trover 
discussion  of  the  uature  of  ponding  of  adsorbates  at  veiy  ,u»  .-.overage  it-)  —  0>  10 
covalent  or  metallic. 
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Ab-Initio  Calculations  of  Seif  Diffusion  Constants  in  Silicon 


P.  E.  Bloechl 

IBM  T.  J.  H'atjon  Research  Laboratory ,  Yorktown  Heights.  USA 

and 

IBM  Research  Laboratory,  Zurich,  Switzerland 


Theoretical  discussion  of  diffusion  mechanisms  has  so  far  focussed  on  ac¬ 
tivation  energies.  Here,  we  report  first-principles  calculations  of  diffusion  con¬ 
stants,  which  allow  a  direct  and  unambiguous  comparison  of  the  relative  con¬ 
tributions  of  two  defect  mediated  mechanisms  for  self-diffusion,  the  vacancy 
mechanism  and  the  interstitial  mechanism.  Each  diffusion  constant  is  the  prod¬ 
uct  of  the  concentration  of  the  native  defect  species  and  its  diffusivity:  we  have 
calculated  the  diffusivity  via  Car-Parrinello  simulations[l]  and  the  concentra¬ 
tion  within  the  local  harmonic  approximation,  that  has  been  suggested  by  LeSar 
et  al.[2'.  The  electronic  contribution  to  the  entropy  of  formation  will  be  dis¬ 
cussed.  Our  results  suggest  that  both  defect-mediated  mechanisms  contribute 
importantly  to  self-diffusion. 

This  work  was  done  in  collaboration  with  W.  Andreoni,  R.  Car,  D.  Laks. 
S.T.  Pantelides,  M.  Parrinello  and  E.  Smargiassi.  It  was  supported  in  part 
by  ONR  Contract  N00014-84-0396,  USARDSG  grant  DAJA  45-89-C0025  and 
CNR  grant  89.00006.69. 

1  R.  Car  and  M.  Parrinello,  Phys.  Rev.  Lett.  55,  2471  (1985). 

2  R.  LeSar,  R.R.  Najafabadi  and  D.  Srolovitz,  Phys.  Rev.  Lett.  63,  624  (1989). 
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Ab-initio  calculation  of  free  energies 
in  real  materials: 

application  to  the  vacancy  formation  in  silicon. 

t.  Smargiassi  and  R.  Car 

International  School  for  Advanced  Studies,  S4014  Trieste,  Italy 

Only  few  ab-imtio  calculations  of  free-energies  of  materials  exist  to  date.  More¬ 
over  these  calculations  have  been  often  based  on  harmonic  or  quasi-harmonic 
approximations.  A  full  inclusion  of  anharmonic  effects  can  be  achieved  by  us¬ 
ing  appropriate  statistical  mechanical  simulation  techniques.  A  popular  method 
consists  in  the  construction  of  a  reversible  path  to  a  state  of  known  free  energy. 
Such  approach  has  been  applied  so  far  only  to  classical  many-bodv  systems  with 
empirical  interatomic  potentials. 

In  this  paper,  we  show  that  this  method  can  be  generalized  to  first-principles 
molecular  dynamics  simulations,  in  which  the  atomic  potential  energy  surface  is 
derived  quantum  mechanically  from  the  electronic  ground-state.  We  apply  our 
scheme  to  compute  the  free  energy  of  formation  of  a  lattice  vacancy  in  silicon  at 
finite  temperature.  This  is  a  basic  quantity  to  evaluate  the  equilibrium  vacancy 
concentration  and  the  diffusion  constants  of  vacancy  mediated  processes.  In  both 
cases  our  results  are  in  good  agreement  with  the  available  experimented  data.  Our 
approach  is  not  limited  to  applications  involving  defective  solids  but  can  also  be 
used  in  a  wider  context,  such  as  the  study  of  phase  transitions  in  materials. 

Support  for  this  work  was  provided  by  the  Italian  Consiglio  Nazionale  delle 
Ricerche  (CNR)  Grant  No.  89.00006.69,  by  the  European  Research  Office  of 
the  U.S.  Army  Grant  No.  DAJA45-89-C-0025,  and  by  the  Italian  Ministry  of 
the  University  and  Scientific  Research  under  the  SISSA-CINECA  collaborative 
project. 
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First  Principles  Molecular  Dynamics  Simulations  of 
Amorphous  Phosphorus  and  Liquid  and  Amorphous  Selenium 

D.  Hohl 

NCSA,  University  of  Illinois  at  Urbana  Champaign 
R.O.  Jones 

IFF,  Forschungszentrum  Julich  (KFA) 

Parameter-free,  first  principles  molecular  dynamics  studies1  (MD/DF)  of  amorphous 
phosphorus  and  semiconducting  liquid  and  amorphous  selenium  are  presented,  focussing 
primarily  on  the  atomic  microstructure  in  these  materials. 

Although  easy  to  prepare  in  various  chemical  and  thermal  ways  and  important  as  industrial 
material,  amorphous  phosphorus  is  poorly  characterized  experimentally  and  theoratically.  A 
simple  cubic  unit  cell  of  constant  volume  comprising  124  P  atoms  was  chosen  to  represent 
the  system  under  consideration,  making  it  the  largest  system  investigated  with  the  MD/DF 
method  to  date.  Starting  from  a  hypothetical  sixfold  coordinated  P,  the  system  is  melted  and 
annealed  to  room  temperature  in  two  cycles.  Trajectories  several  ps  long  serve  to  calculate 
pair  correlation  functions  g(r),  scattering  functions  S(Q),  phonon  densities  of  states  and 
other  static  and  dynamic  quantities.  The  results  serve  to  guide  the  interpretation  of  recent 
detailed  neutron  scattering  experiments  and  are  put  into  perspective  with  our  MD/DF  study 
on  isolated  P„  clusters  in  the  vacuum.2 

The  MD/DF  simulation  of  liquid  and  amorphous  selenium3  in  the  temperature  range 
350-720  K  used  64  atoms  in  a  constant  volume  simple  cubic  cell.  Ring  statistics,  pair 
correlation  functions  g(r),  structure  factors  S(Q),  bond  angle  and  dihedral  angle  distribution 
functions,  phonon  densities  of  states  and  diffusion  constants  D  were  computed  and  compared 
with  experimental  and  theoretical  data  where  available.  Several  structural  models  proposed 
previously  for  Se  are  evaluated  and  their  predictions  compared  with  the  results  of  the  present 
simulation.  Special  emphasis  is  given  to  the  concentration,  structure  and  stability  of  the 
dominant  defects,  and  a  new  structure  and  defect  model  is  presented. 

1  R.  Car  and  M.  Parrinello,  Phys.  Rev.  Lett.  55,  2471  (1985) 

2  R.O.  Jones  and  D.  Hohl,  J.  Chem.  Phys.  92,  6710  (1990) 

3  D.  Hohl  and  R.O.  Jones,  Phys.  Rev.  B  43,  xxxx  (1991) 
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Ab  Initio  Molecular  Dynamics  Simulation  of  liquid  Sodium  - 

Tin  alloy 
by 

G.  Se2  fert 

Institute  tor  Theoretics!  Ph yeics,  TU  Dresden,  Breeder. , 
German  y 

and 

G.  Past or e ,  R .  Car 

S1SSA  -  International  School  for  Advanced  Studies.,  Tries 
Italy 

Abstract 

A  first  principle  mol  ecul  ar-dynami  cs  study  of  liquid  N!sc-  .-,1 
is  presented.  The  structural  properties  (static  structure  fact 
pair  correlation,  bond  angle  distribution)  are  discussed, 
calculated  static  structure  factor  is  in  good  agreement  ». 
that  from  neutron  diffraction  experiments.  A  prepeal  in 
structure  factor  can  be  attributed  to  the  existence  of  tr 
hedral  Sn^  units  in  the  liquid  phase.  The  Sr.fl  com,  1  exes  '  7. 
anions  [Sn^4  J-  are  typical  structural  units  in  solid  h. 
However,  the  calculations  show  also  a  tendency  of  network  fc 
ation  of  three-  and  fourfold  coordinated  tin  atoms  in  the  1 l qu 
As  a  consequence  of  this  network  formation  there  is  an  indi 
1 1  on  of  a  "bi  modal  "  bond  angle  distribution  for  the  tin  at  c>n.  . 
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Alloy  Formation  in  the  Effective  Medium  Theory 
Karsten  W.  Jacobsen 

Laboratory  of  Applied  Physics 
Technical  University  of  Denmark 
DK-2800  Lyngby 
Denmark 

The  formation  of  alloys  are  discussed  within  the  framework  of  the  effective 
medium  theory  for  interatomic  interactions.  The  total  energy  of  simple  metal 
or  noble  metal  alloys  are  shown  to  consist  of  an  electron  density  dependent 
term  which  determines  the  volume  of  the  alloys  and  an  electrostatic  correction 
to  the  atomic  sphere  approximation.  Structural  energy  changes  are  to  a  large 
extend  determined  by  the  latter  term.  The  stabilities  of  the  three  binary 
alloys  formed  by  Cu,  Ag,  and  Au  are  discussed. 

For  transition  metal  alloys  in  the  center  of  the  periodic  table  an  additional 
contribution  to  the  heat  of  formation  comes  from  the  formation  of  an  alloy 
d-band.  This  contribution  is  modelled  using  a  moment  method. 

The  effective  medium  approach  is  shown  to  be  simple  enough  to  allow  for 
analytic  calculation  of  cluster  interactions  from  properties  of  the  constituents 
alone.  Furthermore,  the  method  can  be  used  for  Monte  Carlo  simulations  of 
thermal  properties  of  alloys. 
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Structure  and  thermodynamics  of  Si7Ge,_,  alloys 
from  ab-initio  Monte  Carlo  simulations 

Stefano  de  GironcoIi^Paolo  Giannozzi,1  and  Stefano  Baroni2 

iIRRMA  -  Institut  Romand  de  Recherche  Numerique  en  Physique  des  Matenaux 
PHB  -  Ecublens,  CH-1015,  Lausanne  -  Switzerland 

2SISSA  -  Scuola  Intemazionale  Superiort  di  Studi  Avanzati 
Strada  Costiera  11,  I-SJOIJ  Trieste  -  Italy 

We  present  a  new  method  for  studying  from  first  principles  the  finite-temperature 
properties  of  semiconductor  alloys.  The  alloy  is  considered  as  a  perturbation  with  respect 
to  periodic  virtual  crystal,  and  the  relevant  energies  calculated  by  density-functional 
perturbation  theory.  We  first  show  that — up  to  second  order  in  the  perturbation — the 
energy  of  the  alloy  is  equivalent  to  that  of  a  lattice  gas  with  only  two  body  interactions. 
The  interaction  constants  of  this  lattice  gas  are  particular  linear  response  functions  of 
the  virtual  crystal,  which  can  be  determined  from  first  principles.  As  an  application,  we 
consider  the  case  of  Si7Gei_,.  Ignoring  lattice  relaxation,  the  range  of  the  interactions 
is  very  short.  Atomic  relaxation  renormalizes  the  interactions  and  makes  them  rather 
long  range,  propagating  mainly  along  the  bond  chains.  A  comparison  with  traditional 
self-consistent  calculations  for  small  supercells  demonstrates  that  the  acccuracy  of  the 
perturbative  approach  is  in  this  case  of  the  same  order  as  that  state-of-the-art  density- 
functional  calculations.  Once  the  interaction  constants  have  been  calculated,  the  finite- 
temperature  properties  of  the  alloy  can  be  studied  rather  inexpensively  by  MonteCarlo 
simulations  on  the  lattice  gas.  Our  results  show  that  Sir  Gej  _r  is  a  model  random  alloy  with 
a  miscibility  gap  below  =;  170  K.  The  bond  length  distribution  displays  three  well  defined 
peaks  whose  positions  depend  on  composition,  but  not  on  temperature.  The  resulting 
lattice  parameter  follows  very  closely  Vegard’s  law. 
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Ab-Initio  Lattice-Dynamics  in  Semiconductors 

P.  Pavone0,  S.  Baroni®,  S.  de  Gironcoli6,  and  P.  Giannotri* 

•  Seuola  Internationale  Superiore  di  Studi  Avamati  (SISSA) 

Strada  Costiera  11,  1-34014  Trieste,  Italy 
b  InsUtut  Romand  de  Recherche  Numenque  en  Physique  des  Materiauz  (IRRMA 
PHB-Ecublens,  CH-101S  Lausanne, Switzerland 

We  present  some  calculation*  of  the  lattice-dynamical  properties  of  semiconductors,  based  on  density- 
functional  linear-response  theory.  As  a  first  application,  we  compute  complete  phonon  dispersions  for  the 
elemental  semiconductors  Si  and  Ge,  and  for  the  III- V  semiconductor  compounds  GaAs,  AlAs,  GaSb,  and 
AlSb.  Our  results  are  in  excellent  agreement  with  experiments,  where  available,  and  provide  predictions 
where  they  are  not. 

As  a  byproduct,  we  obtain  real-space  interatomic  force  constants  for  these  materials,  which  are  useful 
both  for  interpolating  the  dynamical  matrices  throughout  the  Brillouin  tone,  and  as  ingredients  of  ap¬ 
proximate  calculations  for  mixed  systems  such  as  alloys  and  mic restructures.  We  study  the  transferability 
of  force  constants  among  ieoclectronic  compounds.  We  find  that  force  constants  of  III- V  compounds  dif¬ 
fering  only  by  their  cation  are  quite  similar,  whereas  this  is  not  true  when  only  the  anion  is  different. 
The  situation  is  intermediate  in  the  case  of  elemental  semiconductors. 

As  a  further  application  we  calculate  the  thermal  expansion  of  Si,  Ge.  AlAs,  and  GaAs  in  the 
quasiharmonic  approximation.  As  intermediate  results  we  obtain  the  Gruneisen  parameters  of  these 
systems  troughout  the  Brillouin  tone.  Also  in  this  case,  the  agreement  with  available  experiments  is 
good. 
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Phonons  and  effective  charges  in  Boron  Phosphide 

H.  W.  L.  Alves  and  K.  Kune 

Laboratoire  de  Physique  des  Solides  associe  au  CNRS 
Universite  P.  et  M.  Curie  4,  PI.  Jussieu  Tour  13  75005  Paris  France 

Static  and  dynamical  properties  of  Boron  Phosphide  are  studied 
within  the  local  density  approximation,  using  the  norm-conservmg 
pseudopotentials  and  plane-waves  expansions.  Calculated  lattice  constant, 
bulk  modulus,  selected  phonon  frequencies  and  Gruneisen  parameters  agree 
with  experiment  to  within  1%  ,5%  for  ag,  Bq  respectively,  and  within  5-7% 

for  the  phonon  frequencies,  whenever  experiment  is  available. 

Real  space  distribution  of  electronic  charge  density  shows  that 
the  valence  charge  is  nearly  equally  divided  between  the  cation  and  anion. 
Previous  estimations  of  this  partition  led,  in  the  past,  to  speculations 
that  the  effective  charges  might  have  the  signs  reversed,  compared  to 
other  lll-V  compounds.  Our  calculations  on  the  quadrupled  suDercelis  wife 
the  longitudinal  displacements  did  not  confirm  the  above  conjecture-at 
the  zero  pressure,  at  least.  Results  corresponding  to  higher  compressions 
will  be  discussed. 


SOFT  PHONONS  IN  CUBIC  BaBi03 

R.  Zeyher  and  K.  Kune 

Max-Planck-lnstitut  fur  Festkdrperforschung,  7000  Stuttgart.  Germary, 
C.N.R.S.  and  University  P.  et  M.  Curie,  75230  Paris ,  France 


The  occurence  of  superconductivity  in  BaPb^Bi^Oj  and  Ba^K^BiOj  has  attractec  ict 
of  interest:  These  systems  share  many  characteristic  features  of  the  Cu-based.  high-7 
superconductors,  but  they  are  considerably  simpler.  The  crystal  structure  of  the  undooed 


BaBi03  derives  from  the  cubic  perovskite  structure  where  the  charge  density  at  the  lattice 
sites  has  been  modulated,  leading  to  a  doubling  of  the  primitive  cell  along  one  axis.  It  has 


been  widely  believed  that  the  distortion  of  the  cubic  perovskite  cell  is  due  to  a  local 
electronic  instability  leading  to  a  charge  disproportion  2Bi4*  -*  Bi3"'  +  Bi5’  with  two 


inequivalent  Bi  sites.  The  aim  of  this  work  was  to  investigate  possible  electronic  and  lattice 


instabilities  of  a  hypothetical  cubic  perovskite  structure  of  BaBi03  and  to  determine 
resulting  distortions  by  means  of  total  energy  calculations. 


We  have  carried  out  the  calculations  by  applying  tne  LDF  theory  in  plane-wave  basis 


and  using  the  norm-conserving  pseudopotentials.  We  verified  that  our  energy  eigenvalues 
tj(k)  correctly  reproduce  the  “band  structures'  obtained  previously  by  other  methods 
(LAPW),  and  checked  that  our  prediction  for  the  lattice  constant  a0  agrees  with  experiment 
(1  %  error).  Keeping  the  ions  at  their  ideal  lattice  positions  we  did  not  find  any  evidence  for 
a  local  electronic  instability  in  cubic  BaBiOj  towards  a  modulated  charge  density  which  was 
observed  experimentally.  Performing  calculations  with  frozen  phonons  at  the  R-point  of  the 
B.Z.  we  found  the  breathing  mode  r2'  to  be  stable.  On  the  contrary,  the  (l.-l.O)-component 
of  the  rotational  mode  r25  is  unstable  in  the  cubic  structure  and  drives  the  transition. 
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a-Gallium:  A  metallic  molecular  crystal 

X.  G.  Gong0,1’,  Guido  L.  Chiarotti”  c,  M.  Pamnello11  d,  E.  Tosatti0  b 

a)  Internationa!  School  for  Advanced  Studies  (SISSA),  Strada  Costtera  11,  1-34014 
Trieste,  Italy 

b)  International  Centre  for  Theoretical  Physics  (ICTP).  P.0.  Box  586,  1-54014 
Trieste,  Italy 

c)  Laboratono  Tecnologie  Avanzate  Superfici  e  Catalisi(TASC),  del  Consorzio  In- 
teruniversitano  Nazionale  per  la  Fisica  della  Materia  (ISFN),  Padnciano  99, 
I-SiOlt  Trieste,  Italy 

d)  IBM  Research  Division,  Zurich  Forschunglaboratorium,  CH-8803  Ruschlikon , 
Switzerland 


Abstract 

We  present  a  first  principle  theoretical  study  of  the  atomic  and  electronic  struc¬ 
tures,  and  of  the  zero  temperature  phases  of  gallium.  The  picture  of  a-Ga  which 
emerges  is  that  of  a  metallic  molecular  crystal,  with  a  strong  Ga2  covalent  bond, 
and  weaker  intermolecular  forces.  It  is  supported  in  detail  by  the  charge  density, 
the  electronic  structure,  and  differential  bonds  stretching  under  pressure.  Anoma¬ 
lous  features  of  a-Ga,  such  as  Knight  shift,  anisotropic  Fermi  Surface  effects,  and 
optical  absorption  find  a  consistent  explanation.  Accurate  X-ray  measurements 
should  reveal  the  Gaj  covalent  bonds. 
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Structural  and  Electron*-  Properties  of  a  —  Sn.  C'dTe  and 
Their  Monolayer  SuperlatticesT 

A  Continenza"  and  A  J.  Freeman 
Department  of  Physics  and  Materiais  Research  Center 
Northwestern  University.  Evanston,  U  60208 

November  8,  1990 


Interest  in  small  direct  gap  semiconductors  is  driving  the  atten':on  on  compounds  based 
on  group  IV  materials  In  particular,  recent  experiments  have  demonstrated  that  the  o-Sn 
zero  band  gap  can  be  opened  by  quantum  confinement  effects  in  ultra-thin  layers  of  o  - 
Sn  deposited  on  suitable  substrates  (such  as  C'dTe).  Moreover,  the  successful  growth  of 
superlattices  of  o-Sn  CdTe  along  different  growth  directions,  has  been  reported  by  several 
groups. 

In  view  of  all  this  interest,  w*  pe-formed  a  systematic  study  of  the  structural  and  ele:- 
tromc  properties  of  Q-Sn.  CdTe  and  their  '001'  monolayer  superlattices  (namely:  SnTeC  d: . 
SnCdTe;.  and  SnjC'dTe)  using  the  total  energy  local  density  all-electron  full  potential  br.- 
•arized  augmented  plane  wave  (FLAPW)  method. 

Our  calculations  reproduce  with  good  accuracy  the  ground  state  properties  of  the  pure 
constituents  and  give  some  insights  regarding  the  properties  of  the  ultra-thin  superlattices 
In  particular,  we  found  that  (  j  the  “compensa'ed"  structure  (i.e  SnjC'dTe)  is  unstable 
with  respect  to  phase  separation,  (ii)  that  the  structure  in  which  Sn  substitutes  for  C'd 
is  the  most  unstable  and  (iii)  tha"  large  tetragonal  distortions  are  caused  by  this  same 
substitution.  In  addition,  in  all  the  structures  studied  the  direct  band  gap  is  favored  over 
the  indirect  T  -  L  gap  and  that  this  direct  band  gap  is  higher  (within  LDA)  than  ir.  pure 
o  -  Sn 


Work  supported  by  the  NSF  (through  the  NT  MRC). 

Present  Address:  Unjversita  degli  Studi  di  L'Aquila  -  Coppito,  1-67010  (Italy). 
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LCAO  Band  Calculation  for  the  Hexagonal  YgFe-iy  and  VgFe-jyN 

Y.  P.  Li,  D.  L.  Weaire,  Hong-Shuo  Li  and  J.  M.  D.  Coey 

Department  of  Pure  and  Applied  Physics,  Trinity  College,  Dublin  2. 
Ireland 

A-bitrasi 

First  principle  band  calculation,  by  using  Orthogonal  Linear 
Conbination  of  Atomic  Orbitals  (OLCAO)  approximation,  have  been 
carried  out  on  Y2Fen  and  the  newly  discovered  Y2Fe  17N  compounds. 
This  iron-rich  rare-earth  nitride  phase  is  expected  to  be  the  main  phase 
for  high-performance  magnets  with  Curie  temperature  and  the 
saturation  magnetization  higher  than  those  of  Nd2Fei4B  magnets.  The 
atomiclike  valence  basis  was  adjusted  according  to  the  local 
environment  of  each  crystallographic  site.  It  is  found  that  the  average 
iron  moment  are  2.40pg  for  YiFep  and  2.50pg  for  YjFenN,  whLh  is  in 

good  agreement  with  Mossbauer  results. 


All  the  correspondences  should  be  addressed  to 
Prof.  Y.  P.  Li 

Department  of  Pure  and  Applied  Physics, 
Trinity  College,  Dublin  2.  Ireland. 

Fax:  +353-1-711759 
Tel:  +353-1-772941 


Relativistic  Atomic  Evaluation  for  the  Exchange 
Interactions  in  the  Rare  Earth  -  Iron  Intermetallies 

Hong-Shuo  Li,  Y.  P.  Li  and  J.  M.  D.  Coey 

Department  of  Pure  and  Applied  Physics,  Trinity  College,  Dublin  2, 

Ireland 

Abstract 

Intra-atomic  exchange  integrals  and  J 4f-6p<  for  rare‘ 

earth  trivalent  ions  were  obtained  from  relativistic  self-consistent  field 
atomic  calculations,  using  the  local-density  approximation  for  exchange 
and  correlation  potential.  The  resulting  effective  (molecular  field) 
exchange  coefficients  nRT  and  nRT  show  a  decrease  in  magnitude  across 
the  rare-earth  series,  which  agree  well  with  those  observed  in  many 
rera-earth  transition-metal  intermetallic  series.  It  is  shown  that  the  4f- 
5d  exchange  interaction  is  far  more  important  than  the  4f-6s  and  4f-6p 
interactions,  as  expected  from  Campbell  model  for  the  exchange 
interactions  occuring  in  rare-earth  transition-metal  intermetallies. 


All  the  correspondences  should  be  addressed  to 
Prof.  Y.  P.  Li 

Department  of  Pure  and  Applied  Physics, 
Trinity  College,  Dublin  2,  Ireland. 

Fax:  +353-1-711759 
Tel:  +353-1-772941 
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Change  of  sign  of  orbital  moment  in  bcc  Fe. 

L.  Severin,  L.  Nordstrom  and  B.  Johansson 

Condensed  Matter  Theory  Group,  Institute  of  Physics,  University  of  Uppsala 
BOX  530  S-751  21  Uppsaia.  Sweden 

We  report  on  electronic  bandstructure  calculations,  including  spin-orbit  coupling,  for 
the  magnedc  transition  metals  Fe,  Co  and  Ni  where  the  spin  moment  is  constrained  to  a 
given  value.  For  bcc  Fe  a  change  of  sign  of  the  orbital  moment  with  decreasing  spin 
moment  is  predicted,  in  contradiction  with  intuitive  thinking.  This  effect  is  not  present 
in  hep  Co  and  fee  Ni.  This  is  explained  in  terms  of  the  detailed  structure  of  the  mj  - 
projected  density  of  states  for  the  bcc  structure. 
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Real-Space  Iterative  Methods  in  Electronic  Structure 
Calculations* 

D.  J.  Sullivan  an|  J.  Bctnholc 

Department  of  Physics,  NC  State  University,  Raleigh,  NC  27695-8202 

Iterative  methods  such  as  the  Car-Parrinello  tfcethod  have  very  successfully  used  plane 
wave  bases  in  the  study  of  bulk  crystals,  surfaces,  liquids,  and  clusters.  However,  since 
plane  waves  provide  uniform  representation  of  fchysical  quantities  throughout  the  unit  cell, 
the  addition  of  a  strong  impurity  requires  al  improvement  of  the  basis  everywhere. 
Consequently,  the  number  of  plane  waves  and  me  cost  of  the  calculation  increase  dramati¬ 
cally  with  the  addition  of,  eg..  a  transirion-me«i  impurity.  In  surface  and  cluster  calcula¬ 
tions,  much  of  the  supercell's  volume  is  outside  the  ionic  core  and  bonding  regions; 
nonetheless,  the  plane  wave  basis  covers  the  vacuum  region  as  accurately  as  the  interior 
region.  Non-uniform  real  space  grids  allow  for  a  local  improvement  of  the  representation  of 
wavefunctions.  i.e.,  in  the  ionic  core  and  bonding  regions,  while  using  fewer  basis  states 
in  the  regions  where  the  wavefunction  is  small  cr  slowly  varying.  We  illustrate  a  real-space 
iterative  method  with  applications  to  atoms  and  diatomics  without  the  use  of  periodic 
boundary  conditions,  and  discuss  the  tailoring  4f  the  basis  and  the  smoothness  of  the  rep¬ 
resentation  fhtr  results  indicate  that  the  prccpdure  is  similar  ir.  efficiency  :c  die  Car- 
Pamnello  method. 

•Supported  hy  ONR.  grant  number  NXCO 1 4-  3  3 1 K  -  04  9 1 . 
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An  LAPW  Approach  to  Separable  First-Principles 
Pseudopotentials 

S.  Goedecker,  K.  Maschke 
IPA,  Ecole  Polytechnique  Federal  e  de  Lausanne 
CH-1015  Lausanne,  Switzerland 

Electronic  structure  calculations  based  on  separable  pseudopotentials  require 
much  less  computer  time  than  calculations  based  on  non-separable  forms.  Un¬ 
fortunately  the  most  reliable  pseudopotentials  sire  non-separable  and  it  is  not 
always  possible  to  make  them  separable  without  loosing  precision.  Following 
ideas  borrowed  from  the  LAPW  method,  we  derive  a  First-Principles  pseu¬ 
dopotential  which  is  separable  by  construction  We  show  that  for  sufficiently 
small  muffin-tin  spheres  it  is  not  necessary  to  recalculate  the  selfconsistent 
potential  in  the  muffin-tin  sphere  at  each  step  of  the  selfconsistency  cycle  as 
it  is  done  in  the  LAPW  method.  This  means,  that  one  reference  potential  in 
the  muffin-tin  sphere  is  enough  to  ensure  good  transferability. 
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A  fast  full-potential  LCAO-method  for  total-energy 
minimizations 

P.  M.  Oppeneer,  W.  Hierse  and  J.  Kfibler 

Institut  fur  Festkorperphysik,  Terhnische  Hochschule  Darmstadt,  D-6100  Darmstadt. 

Federal  Republic  Germany 

A  new  bandstructure  method  has  been  developed,  specifically  for  the  purpose  of  fast 
total-energy  minimizations  of  large  systems.  Most  characteristic  of  this  Linear-Combi- 
nation-of-Atomic-Orbitals-method  is  that  nowhere  position  dependent  meshes  for  nu¬ 
merical  integrations  in  real  space  are  needed.  All  appearing  integrations  are  done 
analytically,  which  makes  the  setting-up  of  the  secular  matrix  fast.  Furthermore,  with 
this  LCAO-method  a  number  of  new  computational  techniques  are  introduced,  namely: 
(a)  an  efficient  and  simple  solution  of  Poisson’s  and  Schrodinger's  equation  for  fuil- 
potentials,  (b)  a  new  treatment  of  the  non-linearity  of  the  LDA  exchange- correlation 
energy,  and  (c)  an  approximation  technique  for  3-centre  integrals.  The  method  de¬ 
scribed  here  is  particularly  suited  for  molecular  dynamics  simulations  and  we  also 
developed  a  scheme  for  this  LCAO-method.  in  which  a  minimum  of  the  total-energy 
is  found  by  iterating  the  electronic  and/or  ionic  degrees  of  freedom. 
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Spurious  pressure  effects  on  LDA 
band  structure  of  semiconductors 

Vincenzo  Fiorentini 

Dipartimento  di  Fisica  Teorica,  Universita  di  Trieste,  Italy 


A  consistent  description  of  the  band  structure  of  semiconductors  in  the  LDA  by  the  pseudopo¬ 
tential  method  (indeed  by  any  scheme  not  reproducing  exactly  the  experimental  equilibrium  vol¬ 
ume)  can  only  be  obtained  by  calculations  done  at  the  theoretical  equilibrium  volume,  0  =  fits- 
The  choice  of  different  cell  volumes  (such  as  the  experimental  volume  Cl„)  is  equivalent  to  a 
spurious  pressure  applied  to  the  crystal. 

These  statements  are  substantiated  by  converged  total  energy  studies,  and  calculations  of  elec¬ 
tronic  states,  pressure  derivatives,  and  effective  masses  of  Si,  Ge,  GaAs,  AlAs  (using  Bachelet- 
Hamann- Schliiter  ab  initio  norm-conserving  pseudopotentials,  Ceperley-Alder  LDA,  and  a  plane 
wave  basis)  :  at  !},/,,  the  discrepancies  in  T  and  L  energy  gaps  are  definitely  smaller  than  usually 
believed;  in  particular,  the  correct  ordering  of  conduction  band  is  obtained  in  Ge.  This  is  easily 
understood  in  terms  of  pressure  derivatives,  whose  calculated  values  agree  well  with  existing 
experimental  data.  Effective  masses  are  also  very  sensitive  to  variations  of  fl,  and,  contrary  to 
common  wisdom,  the  calculated  conduction  effective  masses  in  GaAs  and  AlAs  are  found  to  be 
in  remarkably  good  agreement  with  experiment 
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HARTRKE- FOCK  PLUS  DENSITY  FUNCTIONAL  CALCULATION  OF  GROUND  STATE 
PROPERTIES  OF  SKMI CONDUCTOR  AND  IONIC  COMPOUNDS. 


Mauro  Causa' 

Department  of  Inorganic.  Physical  and  Material  Chemistry. 
University  of  Torino,  via  Gloria  5,  1-10125  Torino.  Italy. 


The  binding  energy  (R.E.).  the  equilibrium  lattice  parameter 
(a0),  the  hulk  modulus  (B)  and  the  centra )- zone  transverse 
optlcnl  phonon  frequency  (v^g!  of  seventeen  diamond  and 
zinc-biend  type  structure  semiconductors  and  of  lithium,  sodium 
and  potassium  oxides  have  been  evaluated.  The  periodic  mb  initio 
Hartree-Fock  I.CAO  program  CRYSTAL  hex  been  used.  Core 
paeudopotentials  have  been  adopted  In  order  to  limit  the 
calculation  to  valence  electrons.  Two  gaueslan  type  function  per 
each  atomic  orbital  plus  polarization  d  function  have  been 
adopted  me  basis  set.  The  quality  of  the  pseudopotcntlal  results 
has  been  checked  by  comparison  with  all  electron  calculations 
performed  on  six  light-atom  systems  (diamond,  silicon.  BN,  BP. 
AfP.  SIC).  The  mean  error  of  the  pscudnpotentia  1  calculations 
with  respect  to  experimental  data  is  -38V.  « 1 . 2V ,  -  G  .  5  V  .  -7.8V 
for  A.K.,  aD .  B  and  v-pg,  respectively.  A  cor  re  1  a t ion- on  1  y 
density-functional  a  posteriori  correction  to  the  Hartree-Fock 
total  energy  1s  performed,  which  reduces  the  nbsolutc  B.E.  mean 
error  to  2.6V. 
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The  effect  of  spin-orbit  coupling  on  ground-state 
properties  of  materials. 


X.  Gonze.  P.C.P.M.,  U.C.L,  Place  Croix  du  Sud,  1,  B-1348  Louvain-la-Neuve,  Belgium, 
and  (present  address)  LASSP,  Cornell  University,  Ithaca.  N.Y.  14853  USA. 

Using  perturbed  density  functional  theory  (PDFT),  we  show  that  the 
modifications  of  ground  state  total  energy  and  charge  density  of  a 
non-magnetic  crystal,  due  to  the  spin-orbit  term,  are  only  second- 
order  effects  in  the  spin-orbit  coupling  constant,  while  this  is  not 
true  for  atoms  or  molecules. 

Then,  we  adress  the  magnitude  of  these  effects  in  As,  Sb,  and  Bi 
(mainly  same  chemical  properties,  but  increasing  atomic  number, 
and  increasing  spin-orbit  coupling  constant).  We  also  try  to  answer 
more  generally  the  following  practical  question  :  given  a  ground 
state  property  to  be  calculated,  what  is  the  magnitude  of  the  error 
caused  by  ignoring  the  spin-orbit  coupling  terms  in  the  calculation  ? 
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First-principles  study  of  graphite  electronic 

properties 

J.-C.  Charlier,  X.  Gonze,  J.-P.  Michenaud. 

University  Catholique  de  Louvain. 

Unite  de  Physico-Chimie  et  de  Physique  des  Materiaux  (PCPM). 
Batiment  Boltzmann,  Place  Croix  du  Sud,  1 
B  -  1348  louvain-La-Neuve  (Belgium). 

The  electronic  properties  of  hexagonal  graphite  have  been  studied  in 
the  framework  of  the  density  functional  technique,  using  non-local 
ionic  pseudopotentials  and  a  large  number  of  plane  waves.  The  band 
structure,  at  the  Fermi  level,  is  compared  with  the 
Slonczewski-Weiss-McClure  model  for  a  three  dimensional  graphite. 
The  valence  charge  density,  the  density  of  states  and  the  Fermi 
surface  are  also  presented. 


X.  Gonze,  P.C.P.M.,  U.C.L,  Place  Croix  du  Sud,  1,  B-1348  Louvain-la-Neuve,  Belgium 
and  (present  address)  LASSP,  Cornell  University,  Ithaca,  N.Y.  14853  USA. 

The  origin  of  the  DFT  band-gap  problem  is  well  understood  :  the 
exchange-correlation  potential  can  change  abruptly  if  one  conduction 
electron  is  added  to  the  solid.  But,  until  now,  the  actual  size  of  this 
effect  has  been  examined  mainly  in  the  restricted  class  of  diamond 
and  zinc-blende  semiconductors,  or  on  model  semiconductors. 

In  the  present  communication,  it  is  shown  that  the  magnitude  of  the 
band-gap  problem  is  strongly  structure-dependent.  From  simple 
theoretical  arguments,  we  explain  why  the  IV-VI  semiconductors 
band  gap  problem  is  weaker,  and  we  relate  this  fact  to  the  good 
agreement  observed  for  the  elemental  semimetal  overlaps, 
especially  Bismuth.  In  the  same  framework,  we  show  that  diamond 
and  Zinc-blende  structures  are  expected  to  generate  large  DFT  band 
gap  discrepancies,  but  small  discrepancies  for  deformation 
potentials. 
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Calculation  of  Gap  Corrections  in  a  Generalized  Density  Functional  Theory 


L.  Fritsche  and  Y.  M.  Gu 

Institut  fur  Theoretische  Physik,  Tl'  Clausthal,  D-3392  Clausthal-Zellerfeld 

As  has  been  shown  by  one  of  the  authors  (L.F.),  the  Kohn  Sham-equations  can  be 
rederived  by  requiring  the  expectation  value  of  the  N-electron  Hamiltonian  to  be  only 
stationary  against  variation  of  the  wave  function,  i.e.  it  must  not  necessarily  attair  a 
global  minimum.  Consequently,  the  entire  KS-scheme  is  valid  not  only  for  the  ground 
state  but  applies  to  excited  N-electron  eigenstates  as  well.  Total  energy  differences  for 
excitations  across  the  gap  of  insulators  and  semiconductors  prove  to  be  expressible  by 
the  respective  band  energy  differences  and  a  correction  term  which  may  be  of  the  same 
order  of  magnitude.  This  correction  term  can  be  approximated  by  a  simple  expression 
which  gives  surprisingly  good  results  for  the  true  band  gaps  in  solid  rare  gases  which 
have  been  chosen  as  a  first  test  case. 
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The  2-D  Hubbard-Peierls  model  In  the  self-interaction  corrected 
density  functional  formalism 

].  A.  Majewski  and  P.  Vogl,  Physics  Department,  Technical  University  of  Munich,  FRG 

We  present  results  for  correlation  functions,  the  energy  gap,  and  the  phonon-induced 
electronic  energy  gain  in  the  2-D  Hubbard-Peierls  model  near  half-filling.  N’oticeab.y, 
we  find  that  lattice  distortions  destabilize  the  antiferromagneuc  ground  state  in  2D 
and  cause  the  electronic  energy  gap  to  decrease.  Contrary  to  the  ID  case,  the  2D 
Hubbard-Peierls  model  thus  exhibits  a  negative  U-tvpe  behavior.  We  have  employed 
both  exact  diagoralizahons  and  the  self-interaction-corrected  local  spin-density  func¬ 
tional  theory  (SIC-LSD).  In  the  latter  case,  the  off-diagcr.al  Langrar.ge  multipliers  are 
taken  into  account  self-consistentlv.  SIC-LSD  gives  excellent  agreement  with  exact  re¬ 
sults  or.  small  lattices;  detailed  comparisons  of  the  exact  and  SIC-LSD  exchange- 
correlatior.  energy  are  presented. 

We  discuss  an  efficient  method  for  minimizing  the  SIC  total  energy  for  a  given  band 
filling  and  spin  configuration.  Various  numerical  mirimizaucr.  techniques  are  com¬ 
pared  to  one  another. 
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OVERVIEW  OF  THE  DEGENERACY-DEPENDENT  SELF-INTERACTION 
CORRECTION  (D-SIC) 

P.  Cortona.  Dipartlmento  dl  Flslca.  University  dl  Genova.  1-16146  Genova.  Italy. 

D-SIC  is  a  method  of  correcting  the  local-density  approximation  to  the  density-functional 
theory  for  the  self-.,  ‘eractlon  effects  which  is  strictly  based  on  the  homogeneous  gas 
theory.  In  this  communication  we  review  the  theoretical  foundations  of  the  method  as  well 
as  his  more  relevant  results.  Particular  attention  will  be  devoted  to  the  comparison  with  a 
variety  of  other  approximations  In  use  in  the  density-functional  theory. 
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Optical  Matrix  Elements  for  Non-local  Pseudopotential  Calculations 


A.  J.  RE  .AD  and  R.  J.  NEEDS 
Cavendish  Lab..  Cambridge  L’mv..  U  K. 

Calculated  absorption  and  emission  spectra  are  important  tools  for  understanding  the 
electronic  structure  of  materials.  The  optical  matrix  element  required  (in  the  dipole 
approximation)  is  usually  calculated  using  the  momentum  operator  £.  but  the  position 
operator  r  can  also  be  used.  These  two  operators  give  the  same  answer  only  if  the  potential 
in  the  Hamiltonian  commutes  with  r,  which  is  not  true  for  non-local  potentialsWe  consider 
the  construction  of  standard  non-local  pseudopotennals  and  show  that  the  r  operator  is  tr.e 
correct  one.  If  a  gauge  transformation  is  used  to  change  the  '  perturbation  into  u  r. 
pcrturbanon  we  find  the  pseudopotential  must  also  be  changed,  and  this  change  included  :n 
the  perturbation.  Numerical  results  for  tree  atoms  are  given  which  confirm  this 
interpretation,  and  show  that  the  errors  caused  by  using  the  n  operator  incorrectly  can  be 
comparable  to  the  matrix  element  itself.  Many  bulk  pseudopotential  calculations  involve 
wavefunctions  and  potentials  expanded  on  a  plane  wave  basis  set  with  repeating  boundary 
conditions.  We  present  a  direct  method  for  calculating  the  r  matrix  element  for  these 
systems,  allowing  the  correct  oprcal  matrix  element  to  be"  calculated  using  standard  density 
functional  theory  wavefunctions. 
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A  SEPARABLE  FINITE  GAP  POTENTIAL  METHOD  AS  NE'.V  METHOD 
FOR  CALCULATION  OF  ELECTRON  ENERGY  SPECTRUM  OF  SOLI IS 

Bar’yaktar  V.G.,  Belokolos  E.D.,  Korostil  A.M. 

Institute  for  Physics  of  Metals  of  the  Ukrainian  cademy  of 
Sciences,  Kiev  State  University,  Kiev,  USSR 

A  new  method  for  calculation  of  the  electron  energy  spect¬ 
rum  of  solids  that  based  on  using  of  exact  expression  for  se¬ 
parable  potential  and  wave  functions  of  electron  system  with 
finite  set  gaps  in  spectrum  is  suggested.  The  method  gives  a 
possibility  to  describe  with  simple  formula  and  analyse  in  de¬ 
tails  all  properties  of  energy  electron  spectrum.  In  particu¬ 
lar  methods  of  free  electrons  and  strongly  bound  electrons 
are  shown  to  be  limit  cases  of  the  proposed  method. 

Within  the  framework  of  separable  one  gap  potential  we  have 
calculated  an  electron  energy  spectrum,  an  isoenergetical  sur¬ 
faces,  an  electron-phonon  interaction  constant,  etc.  and  have 
applied  these  results  for  explanation  different  properties  of 
conductors  and  superconductors. 
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HAMILTONIAN  CHAOS  IN  Tuc  ENE^Oy  i  R ANS^Eh  OS  A  NONLINEAR 


0,  HENNIG  are  B.  ESSER 

Facnoereicn  Physik  aer  numboiot-Universitat  zu  Berlin,  Ins 
Theoretische  Physik 
Invaiidenstr.  42,  0-1040,  Berlin.  FRG 
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Electronic  Structure  of  Conducting  Copolyxers 

Nagendra  N.  Beladakere,  R.  S.  Rao  and  E.  D.  Nalhctra 
National  Physical  Laboratory 
Dr.  K.  S.  Krishnan  Road 
New  Delhi.  110012.  INDIA 

Abstract 

There  is  a  quest  for  electrically  conducting  copolymers  and 
blends  as  it  is  possible  to  synthesize  a  single  polymeric 
material  possessing  both  electrical  conductivity  and  desired 
mechanical  properties.  It  is  both  useful  and  important  to 
study  the  electronic  structure,  moieties  and  disorder  of 
these  materials.  The  valence  effective  Hamiltonian  technique 
has  been  used  by  us  to  calculate  the  elecctronic  structure 
of  poly ^naphthalene  oxide  pyrrole  and  polyanisidine  aniline 
copolymers.  Our  calculated  results  are  compared  with  the 
experimental  results.  Attempts  are  being  made  to  study  these 
materials  as  disordered  systems. 
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Some  improvements  for  EMT,  Application  to  O  in  YBa;Cu30r 

Andre  Zwartkruis,  Institute  of  Theoretical  Physics, 

Chalmers  University  of  Technology,  S-41296  Goteborg,  Sweden 

We  have  calculated  diffusion  paths  for  interstitial  0  in  YBajCujO?  by  using  the 
Effective-Medium  Theory  for  chemical  binding  (the  version  described  in  J.K.  Norskov, 
Phys.  Rev.  B26  (1982)  2875).  We  have  found  that  for  internal  consistency  of  the  theory, 
one  needs  to  include  an  extra  kinetic  repulsion  term  which  can  be  derived  in  a  systematic 
way  from  Density  Functional  Theory.  For  larger  impurities,  as  e.g.  oxygen  atoms,  the 
contribution  from  the  extra  term  is  large,  and  unphysical,  deep  wells  in  the  potential 
energy  surface  are  removed.  Inclusion  yields  two  further  advantages.  Firstly,  one  can  now 
easily  establish  the  equivalence  of  EMT  and  another  embedding  scheme,  CEM.  recently 
proposed  by  Kress  and  DePristo  (J.  Chem.  Phys.  87  (1987)  4700).  Secondly,  one  can 
justify  to  much  larger  extent  the  splitup  of  the  first-order  corrections  in  either  the  EMT 
or  the  CEM  scheme  as  a  sum  of  pairpotentiais. 

Relaxation  was  included  in  the  computations.  For  0  diffusion  paths  along  the  c-direc- 
tion  we  find  barriers  of  about  1.5  [eV]  and  for  paths  in  the  b-direction.  parallel  to  the 
Cu-0  chains,  we  find  barriers  of  typically  0.5  [eV], 
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Light  lanthanide  impurities  in  iridium:  from  non-magnetic  4/ states  to 
large  local  magnetic  moments  with  both  a  spin  and  an  orbital  part. 

Lars  Nordstrom  and  Borje  Johansson 
Phys.  Dept.,  Uppsala  Univ.,  Box  530,  Uppsala.Sweden. 

The  electronic  structure  of  light  lanthanide  impurities  in  the  small  volume  host  iridium  is 
calculated  by  the  LMTO  Green's  function  method.  In  the  calculation  both  spin-  anc 

orbital  polarization  are  allowed  to  occur.  Ce  and  Pr  are  found  to  not  form  local  magnetic 
moments  while  Nd  and  Pm  form  local  magnetic  moments  with  both  the  spin  and  the 
orbital  pan  almost  saturated.The  hybridization  between  theiridium  5 d  states  and  the 
impurity  4 f  states  with  t2u  symmetry  is  shown  to  be  the  mechanism  for  the  unstable 
magnetic  moments  in  Ce  and  Pr.  The  results  are  compared  with  and  are  found  to  be  in 
agreement  with  data  from  experiment. 
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*Ab  initio  total-energy  calculations  on  hydrogen  ir.  alimirgM* 

A.  Do  Vito  and  N.J.  Cillan 
Physics  Department,  Unlverstiy  of  Keels, 

Keele,  Staffordshire,  S?5  5BG,  U.K. 


He  have  performed  pseudopoter.tial  total-er.ergy  calculations  on 
hydrogen  in  aluminium  using  the  Kleintna.n-By  lander  form  of  the  Baohcl: 
Haraann-Schiuter  oseudcpctential  for  aluminium  and  the  bare  Couiomfcic 
potential  for  the  proton.  The  physical  quantities  studied  arc  the 
hydrogen  heat  of  solution,  the  relative  energies  of  different  inter¬ 
stitial  sites,  the  diftusivnal  migration  energy  and  the  hydrogen- 
vacancy  binding  ar.eriy.  In  all  cas?s,  the  calculations  fully  include 
the  relaxation  of  the  alur.miur  ions.  We  shall  present  a  comparison 
cr  the  calculated  energ-u  with  available  experimental  results, 
and  mans  of  the  electron  density  distribution  surrounding  the  hydroge 
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Quantum  Behaviour  of  Muons  in  Crystalline  Silicon 
C.  Pennetta 

Dipartimento  di  Fisica,  Universita  di  Lecce,  Via  Arnesano,  I-7S100  Lecce 

A.  Baldereschi 

Dipartimento  di  Fiiica  Teoriea,  Universita  di  Trieste,  Strada  Costiera,  I- 34! 00  Trieste 
IRRMA,  PEB  -  Ecublens,  CH-1015  Lausanne 

We  report  the  results  of  a  fully  quantum-mechanical  calculation  of  the  states  of 
an  isolated  positive  muon  in  crystalline  Si.  The  interaction  of  the  muon  with  the 
host  is  described  with  an  effective  potential  derived  within  the  Born-Oppenheimer 
approximation.  The  one-particle  Schrodinger  equation  of  the  muon  is  solved  using  a 
localized  orbital  basis  set.  We  find  that  muons  do  not  behave  like  protons,  as  commonly 
assumed.  In  particular,  even  in  the  ground  state  the  muon  wave  function  is  significantly 
delocalized.  Furthermore,  the  general  features  of  the  energy  spectrum  and  of  the  excited 
state  wave  functions  indicate  that  quantum  effects  strongly  affect  the  behaviour  of 
muons  inside  the  crystal. 
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EXCITONS  BOUND  TO  AN  IONIZED  DONOR  IMPURITY  IN 
SEMICONDUCTORS 


J.J.De  Groote,  C.D.  Caldwell*,  O.  Hipolito,  M.H.  Degani  and  J.E.M.  Hornos 
Departamento  de  Fisica  e  Ciencia  dos  Materials 
Instituto  de  Fisica  e  Quimica  de  Sao  Carlos 
Universidade  de  Sao  Paulo,  C.P.  369 
13560  -  Sao  Carlos  -  Sao  Paulo  -  Brasil 

’  University  of  Central  Florida 
Orlando,  Florida  32816  -  U.S.A. 


Abstract 


Hyperspherical  coordinates  are  used  to  investigate  the  three-particie  complex 
corresponding  to  an  exciton  bound  to  a  Coulomb  center  in  a  semiconductor.  The  ground 
and  excited-state  potential  curves  and  binding  energies  are  calculated  as  function  of  the 
mass  ratio  of  the  hole  and  electron.  The  results  we  have  obtained  are  in  very  good 
agreement  with  experimental  data  for  several  semiconductor  materials.  The  technology 
is  extend  to  more  complex  systems  as  biexcitons  and  clusters  of  excitons.  The  solution  of 
the  partial  equations  in  this  case  is  impossible  by  tradicional  methods  and  a  procedure 
based  in  the  simulated  annealing  is  sugested. 
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THEORETICAL  INVESTIGATION  OF  THE  MOTIONAL  STATES  OF  THE 
NITRITE  ION  IN  ALKALI  HALIDE  CRYSTALS 

G.  K.  Pandey,  M.  Massey  and  K.  L.  Pandey 
Physics  Department,  Allahabad  University,  Allahabad,  India 

A  microscopic  model  for  the  subsitutional  NO^  impurity  ion 
in  NaCl-structure  alkali  halide  crystals  is  developed  on 
the  basis  of  a  specific  theoretical  potential.  Coulomb, 
induction,  and  short-  range  interactions  are  included  to 
obtain  expressions  for  the  barriers  hindering  the 
rotational  motion  of  the  impurity  in  the  matrix  cage. 
Preferred  o r i en tat i ons ,  off-center  displacements,  and 
barrier  parameters  for  this  impurity  are  calculated  in 
KC1,  KB r  ,  and  KI  host  crystals.  The  results  are  used  to 
explain  most  of  the  features  of  the  infrared  spectra  in 
these  impurity  systems. 

CN.B.  Accepted  for  publication  in  Phys.  Rev.  B  3 
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The  total  energy  pseudopotential  method  used  on  stacking 
faults  in  Aluminium 


B.  HAMMER 

Laboratory  of  Applied  Phyiscs,  Technical  University  of  Denmark 
Building  307,  DK-2800  Lyngby,  Denmark 
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MANY-BODY  EFFECTS  W  A  SCHOTTKY  BARRIER 

J.  Charlesworth,  R.W.  Godby  and  R  J.  Needs 

Cavendish  Laboratory.  University  of  CambridgeMadiagley  Rd.,  Cambridge  CB3  OHE.  UK 

Research  into  the  rectifying  action  of  the  Schottky  barrier  has  been  going  on  since  the  1930' s  but  there  still 
exists  no  accurate  theory  to  explain  the  Schottky  barrier  height  or  the  relative  importance  of  different  effects. 
Recently,  calculations  have  been  done  on  a  range  of  interfaces  using  the  LDA  in  DF7.  but  these  greatly 
underestimate  the  band  gap  in  the  semiconductor.  In  addition,  the  LDA  calculations  omit  important  correlation 
effects,  such  as  the  image-like  band-gap  narrowing  described  below. 

We  investigate  many  body  effects  in  the  Al-GaAs(llO)  Schottky  barrier  using  ab-initio  computational  many 
body  theory  in  the  form  of  the  GW  self  energy  where  the  Greens  function.  G,  and  the  screened  Coulomb 
interaction,  W,  may  both  be  gained  from  LDA  approximation  calculations.  This  method  is  known  to  give 
accurate  band  gaps  for  a  number  of  materials  and  so  allows  an  accurate  evaluation  of  the  barrier  heights. 

We  have  extended  the  usual  treatment  to  allow  us  to  examine  the  variance  in  the  self  energy  contribution  to  the 
many  body  Hamiltonian’  as  a  function  of  distance  from  the  Al-GaAs  interface  which  we  may  consider  to  be  a 
state  dependent  local  potential  which  gives  the  effective  exchange  and  correlation  potential.  Because  the  self 
energy  is  integrated  over  all  space,  this  effective  potential  in  the  semiconductor  includes  the  correlation  of  an 
electron  with  those  in  the  metal.  In  particular,  it  includes  a  image-type  potential  term,  which  being  a  measure 
of  the  election  correlation,  is  absent  in  local  density  calculations. 

By  evaluating  this  effective  potential  for  various  states  over  equivalent  regions  (Le.  over  slabs  whose  length  is 
commensurate  with  that  of  the  semiconductor  cell  parameter),  we  can  evaluate  the  importance  of  the  image-type 
terms  as  a  function  of  the  distance  to  the  interface.  We  thus  observe  the  variation  of  the  effective  potential  and 
so  measure  the  resultant  band  bending  and  band  gup  narrowing  for  both  valence  and  conduction  bands. 
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Pseudopotentia!  calculations  of  the  valence  band  offsets  at 
the  ZnSe/Ge,  ZnSe/GaAs  and  GaAs/Ge  (110)  interfaces: 
effects  of  the  Ga  and  Zn  3d-electrons 

A.  Qteish  and  R.  J.  Needs 

Cavendish  Laboratory,  Madingley  Road,  Cambridge  CBS  OHE,  U.K. 


Valence  band  offsets  at  semiconductor  interfaces  involving  Zn  and  Ga  atoms  are  cal¬ 
culated  using  a  first-principles  pseudopotential  technique.  We  show  that  accurate  re¬ 
sults  can  be  obtained,  with  the  3d-electrons  treated  as  core  states,  provided  that  non¬ 
linear  core  exchange-correlation  corrections  axe  included.  In  this  work  the  ZnSe/Ge, 
ZnSe/GaAs  and  GaAs/Ge  (110)  interfaces  have  been  considered.  The  calculated  va¬ 
lence  band  offsets  are  in  excellent  agreement  with  the  experimental  data  and  with  linear- 
muffin-tin-orbital  calculations  (in  which  the  3d-electrons  are  treated  as  valence  states!. 
Furthermore,  we  show  that  such  effects  can  be  incorporated  within  the  model-solid 
theory  of  Van  de  Walle  and  Martin. 
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INTERFACE-LIMITED  ORDERED  COMPOUND  FORMATION: 
APPLICATION  TO  SI/GE 

K.A.  Madtr1  5,  E.  von  Kind1,  and  A.  Baldereschi 2,3 

1  Lab.  fur  Festkorperphysik,  ETH,  CE-8093  Zurich,  Switzerland 

2  Inst.  dt  Physique  Appliquee,  EPF,  CE-1015  Lausanne,  Switzerland 
3IRRMA,  PE  Ecublens,  CE-1015  Lausanne,  Switzerland 

A  static  model  for  the  atomic  structure  of  semiconductor  interfaces  which  takes  into  ac¬ 
count  both  strain  and  chemical  energies  is  presented.  It  is  shown  that  space-confined 
ordered  intermixing,  limited  by  enthalpy  rather  than  diffusion,  can  occur  at  the  interface 
in  spite  of  bulk-instability  of  the  new  phase.  We  apply  the  model  to  Si/Ge  heterojunctions 
and  strained-layer  supsrlattices,  where  strain  and  chemical  energies  are  of  the  same  order 
of  magnitude.  The  qualitative  trends  of  our  results  agree  with  experimental  data,  but  the 
occurrence  of  ordered  SiGe  phases  at  high  temperature  remains  still  unexplained. 
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VALENCE  BAND  OFFSET  AT  Si/Ge  INTERFACES 

L. Colombo^1  \  R.Resta^2)  and  S.Baroni*2) 

O)  Dipartimenio  di  Fisica,  Umversitd  dt  Milano,  via  Celona  16,  201  S3  Milano 
(?)  SISSA  -  Strada  Costiera  11,  34014  Trieste 

The  Si/Ge  superlattices  (SL's)  are  the  paradigmatic  ca»e  of  isovalent,  non- polar,  lattice  -mismatched 
system.  In  this  work  we  consider  thick  (bn1  still  pteudomorphic)  SL’s  such  that  each  of  the  Si  Ge 
interfaces  can  be  considered  as  isolated,  and  the  main  issue  is  whether  the  valence  band  offset  (VBO  is 
tunable  by  means  of  either  chemistry  (t.e.  altering  the  abruptness  of  the  interface)  or  strain.  The  present 
theoretical  study  is  performed  within  the  state-of-the-art  DFT-LDA  framework. 

As  far  as  chemistry  is  concerned,  our  finding  is  that  the  VBO  is  independent  of  chemical  composition 
and  abruptness  of  the  interface.  Detailed  calculations  for  (OOl)-oriented  heterojunctions  show  that  the 
difference  in  VBO  between  a  sharp  Si/Ge  interface  and  a  mixed  one  is  smaller  than  our  present  compu¬ 
tational  noise  (0.01  eV).  Similar  results  are  expected  for  different  orientations  as  well,  being  the  Si  Ge 
interface  isovalent. 

As  far  as  strain  is  concerned,  our  main  finding  is  that  the  VBO  at  strained  Si^Gt  junctions  :s  not 
affected  by  interface  features  (neither  strain  nor  buckling),  while  it  does  depend  upon  the  macroscopic 
strain  conditions,  governed  by  bc-#h  substrate  and  orientation.  The  nontrivial  orientation  dependence  of 
the  strain-induced  VBO  is  a  general  qualitative  features  of  the  Absolute  Deformation  Potentials  i  ADP 
The  VBO  for  (001)  interfaces,  when  expressed  in  terms  of  the  alloy  concentration  of  the  5irGe;_j 
substrate,  takes  the  form  (in  eV): 

VBO(x)  =  0.42**0.46(1-*), 

and  displays  perfect  linearity.  The  tunabLity  feature  of  the  VBO  is  therefore  present,  although  in  a  very 
limited  range  (0.04  eV),  due  to  the  very  close  values  of  the  ADP’s  for  pure  5i  and  Ge 

jl'  R.Resta,  L. Colombo  and  S.Baioni,  Phys.Rev.  B41,  12358  (1990) 
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TUNING  BAND  OFFSETS  AT  SEMICONDUCTOR  HETEROJUNCTIONS 

BY  INTRALAYER  DEPOSITION 

M.  Peressi,“  6  S.  Baroni,6  R.  Resta,6  and  A.  BaldereschU  r 

(a)  £)ffi  Untversild  di  Trieste,  Strada  Coshera  11,  1-34014  Trieste  -  Italy 

,6)  SISSA,  Strada  Costiera  11,  1-34014  Trieste  -  Italy 

,c)  IRRMA,  PHB  -  Ecublens,  CB-1015,  Lausanne  -  Switzerland 

Band  offsets  can  be  modified  at  temiconductor  hetero junctions  or  created  at  homojunctions 
by  depositing  thin  intralayers  of  atoms  with  different  chemical  valence.  Only  very  recently 
this  has  been  realised  in  the  case  of  Ge/GaAs/Ge  (111),1  and  GaAs/Si/AlAs  (001)2,  and  the 
corresponding  modifications  of  the  band  offsets  have  been  measured. 

We  present  here  a  theoretical  study  (based  on  density-functional  self-consistent  and  linear 
response  techniques)  of  Si  and  Ge  intralayers  deposited  at  GaAs/AlAs  (001)  heterojunctions, 
and  at  GaAs  and  AlAs  (001)  homo  junctions.  Our  results  show  that  the  band  offset  can  be  tuned 
in  a  wide  range  (about  1  eV),  being  very  sensitive  to  the  nominal  thickness  and  the  abruptness 
of  the  intralayer.  A  comparison  with  recent  experiments  for  GaAs/Si/AIAs  suggests  that  in  the 
real  samples  Si  atoms  are  confined  over  two  atomic  planes  for  coverages  lower  than  about  0.5 
monolayers,  whereas  for  higher  coverages  Si  diffusion  occurs. 

[1]  J.T.  McKinley,  Y.  Hwu,  B.E.C.  Koltenbah.  G.  Margaritondo,  S.  Baroni,  and  R.  Resta.  J. 

Vac.  Sci.  Technol.  B,  in  press. 

[2-  t.  Sorba,  G.  Bratina,  A.  Antonini,  J.F.  Walker,  M.  Mikovic,  G.  Ceccone  and  A.  Franciosi. 

Phys.  Rev.  B,  in  press. 
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TUNING  EA$£)  3  AN  I  DISCONTINUITIES  IN  StMICCNCUCTCR  HETEF.GJUNCTI  i: 

A.  Munoz 

Departamento  de  Fisica  Fundamental  y  Experimental,  rac.  Fi s  ica  .  I'mvers  i  cad 
de  La  Laguna,  E38204  LA  LAGUNA,  TENERIFE.  SPAIN 

The  possibility  of  tuning  ^eterojunctions  band  offsets  has  been 
recently^  theoretically  proposed  ,  and  has  been  made  experimentally 
possible  .  he  idea  in  this  studies  is  to  generate  a  local  dipole  at  the 
interface  region,  this  dipole  can  be  different  depending  on  the  interface 
composition  and  geometry,  so  that  provide  a  mechanism  for  the  possibility 
of  band  oftsecs  engineering  by  using  thin  layers  at  interface,  that  suooort 
strongly  the  idea  that  m  some  cases  che  extrinsic  charge  neutrality  level. 

In  this  work  we  will  present  different  results  for  tne  case  of  tmr.  Si 

layers  embedded  in  Ce/CaAs  (100!  and  GaAs/AlAs  (10C;  he  t  ero;  ur.c  t  ions  .  us  mg 
scare  of  the  art  calculations  based  or.  local  density  functional  rherrv. 
using  nonlocal  normconservir.g  pseudopotentials  .  in  this  case  we  car.  extra:: 
not  only  the  banc  ctcsets  ,  we  also  provide  the  format. cr.  entr.altv  fcr  c.nis 
interfaces. 

This  work  was  supported  by  the  Cobier.no  Autonot.o  Car.sn: 

1-A.  Munoz.  N.  Chetty  and  R.M.  Martin  .  Phys .  Rev  B  41  ,  Li 

Flores.  A.  Munoz  and  J.C.  Durar..  Appl.  Surf.  Sci.  41/42  1- 

Perez.  A.  Munoz  and  F.  Flores.  Surf.  Sri.  226  3".  { 1 9 8 c 
'  '  U  Era  t  ir.a  ,  L.  Scrba.  G.  Ceccone,  A .  Antonin  i  .  M .  Micovi: 

F.  Tomcr.as mi  and  A.  Franciosi,  ICth  ICPS  and  EIGS3-  l: 
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OP'ICAu  PROPERTIES  OR  SUPER_A“T  I CE  5 

C.~se*-03-..  t-  . *"  .  Rc  . atcc  cj  .  G.Tnecccrcu 

Aristotle  University  cf  Thessaloniki,  Physios  Department 
54006  Tnessa lohi k i ,  GREECE 


In  recent  years,  the  ootical  orocerties  of  the  supe^Iatt: ces 
( Su«  >  have  received  a  Great  deal  of  attention.  In  this  wc,’<  we 
show  that  the  dielectric  -function,  for  thick  Sl„.  becomes  fe 
average  of  tne  dielectric  functions  o*  the  materials  that  :t 
sists  of.  In  order  to  show  this,  we  decompose  t-e  E^  wave;;-ct i on 
in  two  darts.  one  cart  whicn  has  amolitude  in  the  f-.r  =  t 
stituent  and  another  one  which  has  amdituce  in  the  second  con¬ 
stituent.  By  using  tnis  decomposition  ,  the  dielectric  functic- 
Decomes  a  sum  of  di  “ere-t  terms.  It  is  snown  that  most  o-  tnen 
are  generally  small  ana  only  two  terms  a^e  important  t  -  .  r  - 

Sl_«.  These  terms  coi'-sstor;  to  tne  dielectric  functic-  c -  t-e 
bulk  constituents.  Tr.-.s  decomposition  is  also  useful  for  thin  s_- 
oerlattices.  In  tnis  case  it  is  possible  to  fine  the  origin  o- 
tne  different  structures  which  aopear  in  tne  dielectr; 

As  a  numerical  example,  we  apciv  the  above  decomposi 
GaAs/AIAs  SL. . 


c  'jnctior,. 
t  ion  tc  tne 
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FIRST-PRINCIPLES  CALCULATION  OF  THE  MAGNETOCRYSTAL  I  P--H 
ANISOTROPY  ENERGY  OF  Co.X„,  METALLIC  MULTILAYERS. 

G.H.O.  Daalderop.  P.J.  Kelly  and  M.F  H  Schuurmans. 

Philips  Research  Laboratories.  5600  JA  Eindhoven,  The  Netherlands 

The  results  of  first  principles  calculations  of  the  magnetocrystalline  anisotropy  energy 
(MAE)  of  (11!)  metallic  multilayers  Co^„  (X  =  Cu.  Pd  and  Ag;  n  m  ?  or  6 
are  presented.  The  MAE  calculations  were  carried  out  within  the  local-spin  density 
approximation  using  the  LMTO  method  For  Co„Pd„  the  easy  axis  is  found  to  be 
perpendicular  to  the  superlattice  plane,  in  agreement  with  experiment.  The  anisotropy 
energy  decreases  with  increasing  n  and  is  not  strongly  dependent  on  m\  in  our  calcu¬ 
lations  Co,Pd2  and  Co,Pds  have  approximately  the  same  anisotropy  energy.  The  easy 
axis  is  found  to  be  perpendicular  to  the  superlattice  plane  in  the  series  Co,X:  with  the 
maximum  anisotropy  energy  occuring  fer  X  =  Pd.  By  neglecting  the  Pd  on-sitc 
spin-orbit  coupling,  the  anisotropy  energy  is  reduced  by  a  factor  ~2.  Th;s  underline1 
the  importance  of  having  not  only  low  symmetry  but  ako  large  spi r -orbit  coupling 
By  calculating  the  MAE  for  bulk  trigonal  Co  for  a  large  range  of  r,,j  ratios  \n;  car 
show  that  the  large  MAE  found  above  results  from  the  presence  of  ihc  Co  PJ  inter¬ 
face. 


85 


Electronic  structure  of  metallic  superlattices:  Mo/V 

S.  Papadia 

Institute  of  Theoretical  Physics, 

Chalmers  University  of  Technology,  S-412  96  Goteborg,  Sweden 

K.  Karlsson  and  P.  O.  Nilsson 

Department  of  Physics, 

Chalmers  University  of  Technology,  S-412  96  Goteborg,  Sweden 
Abstract 

Electronic  structure  calculations  have  been  performed  for  a  metallic  superlattice 
system  (Mo/V)  using  the  linear  muffin-tin  orbital  method  (LMTO)  in  the  atomic 
sphere  approximation  (ASA).  Total  energies  have  been  calculated  within  the  density 
functional  formalism,  making  the  local  density  approximation  for  the  exchange- 
correlation  part  of  the  energy  functional.  The  aim  is  to  explore  the  band  structure  as 
well  as  the  geometric  structure  of  such  systems  and  to  determine  the  electron  density 
variation  in  order  to  understand  hydrogen  storage  in  these  materials.  Results  are 
presented  for  superlattices  consisting  of  1,  3,  and  5  monolayers  of  each  constituent. 
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Estimation  of  the  ideal  fracture  strength  for  the  jellium 
model  and  rigorous  theorems  for  bijellic  interfaces 


P.  Ziesche  and  R.  Kaschner 

Technische  Universitat  Dresden,  Sektion  Physik,  Germany 

Abstract 

Recent  considerations  of  the  authors  concerning  the  adhesion 
between  solids,  especially  a  new  interpolation  formula  for 
the  cleavage  force  as  well  as  the  numerical  solution  of  the 
coupled  Euler  and  Poisson  equations  of  the  gradient  expan¬ 
sion  method  applied  to  the  jellium  model  of  cleavage  are 
summarized.  Rigorous  theorems  for  the  bijellic  interface 
are  (i)  used  to  test  the  calculations,  and  (ii)  generalized 
to  the  case  of  finite  slit  width.  The  results  include  data 
for  the  cleavage  force  as  a  function  of  separation  and  den¬ 
sity,  the  linear  and  quadratic  force  constants  (which  are 
related  to  linear  and  quadratic  response  functions),  and  the 
interfacial  energy. 
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Structure  and  Dynamics  of  Laif>e  Clusters:  the  Icosahedral- 
Cuboctahedral  Transition  in  El  Clusters  and  Solid  C$q  * 

Jae-Yel  Yi.  Q.  Zhandand  J.  Bemholc 
Department  of  Physics,  N'ortMCarolina  State  University 
Raleigh,  NC  2t695-8202 


The  13-  and  55-atom  clusters  are  among  the  smaller  clncterc  for  which  both  cubcstehtdrol  zr.i 
icosahcdral  structures  exist.  We  have  investigated  ft  relative  energetics  of  the  rwo  structures  for 
A1  usine  the  Car-Parrinelln  method  tf  only  tfeathirg  mode  relaxation:  are  allowed,  the 

•  •« •  •  •  ••••_  ‘AVViT  w'MwCvtuTivwuv/u  ptCiw-nM  10* 

Using  especially  developed  procedures  for  efficientsimulated  annealing,2  we  have  transformed  the 
ideal  cuboctahedra  and  icosahedra  to  distorted  icosartedra.  The  icosahedra  are  thus  the  ground  state 
structures,  although  the  energetic  difference  between  these  structures  and  "weakly"  annealed 
cuboctahedra  is  relatively  small.  The  calculated  ionization  potentials  and  electron  affinities  are  in 
much  better  agreement  with  experiment  after  annewing,  which  demonstrates  the  importance  of  in¬ 
cluding  structural  relaxations  in  calculations.  For  fie  55-atom  clusters,  the  structure  factors  of  a 
"weaklv"  annealed  but  definiteiv  cuboctahedral  stricture  and  the  fully  annealed  icosahedron  are 
similar 'and  quite  different  from  either  of  the  ideal  fructures.  Therefore,  a  structural  classification 
of  small  clusters  on  the  basis  of  indirect  measuremejts  is  difficult  and  not  necessarily  reliable. 


We  also  report  results  of  Car-Parrinello  simulations  for  the  recently  discovered  Cso  molecular 
solid  2  which  is  »  new  and  unexpected  form  of  soltl  carbon.  Symmetry  unconstrained  opr.rn.aa- 
tion  of  the  structure  results  in  alternating  double  and  single  bonds,  similar  to  the  ones  predicted  for 
the  gas  phase  molecules.  The  individual  bond  angl$  deviate  by  at  most  1°  from  their  ideal  values 
of  108°  and  120.°  The  computed  density  of  states  agrees  very  well  with  very  recent  photoe mission 
results.4  Although  the  diameter  of  each  Qo  fullereje  is  7.2  A,  charge  density  contours  show  that 
the  imermolecular  interactions  are  small  at  the  experimental  lattice  spacing  (10.04  A).  Dynamical 
simulations  reveal  that  each  molecule  rotates  at  terqperamres  as  low  as  -300  K.  At  T=  1500  K  the 
primary  motion  is  still  rotation,  with  the  distortion!  from  the  T=0  structure  being  small.  The  las: 
result  indicates  that  fullerenes  should  be  stable  at  »ery  high  temperatures  and  is  consistent  with 
their  formation  under  verv  hot  conditions  in  laseisahlarinn  and  grapKiw  arc  sour:*'  From  '-he 
T=1500  K  simulation  we  have  computed  the  vibratonal  frequency  spectrum  at  this  temperature. 
The  peaks  in  the  speenrum  are  in  good  agreement  wit;  infrared  absorption  data. 

•Supported  by  ONR,  gram  number  N00014-88-K-4491. 
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ENERGETICS  AND  SEGREGATION  IN  NA/CS  HETEROCLUSTERS 

MJ.  L6pez,  A.  Mertenes'',  M.P.  Irtiguez,  J.A.  Alonso 
Departamento  da  Ffatca  Tedrice  Atdmica,  Molecular  y  Nuclear,  Unlvereldad  de 
Valladolid,  Valladolid,  Spain. 

*  Departamento  da  Flalca  Modema.  Unlvereldad  da  Cantabria,  Santander,  Spain 
The  equilibrium  geometry  and  the  corresponding  ground  state  energy  ere 
calculated  (or  NamCe(l  clustert  with  composition  man  end  m=2n,  by  minimizing 

the  total  cluster  energy  with  respect  to  both,  ionic  positions  end  the 
electronic  density.  Density  functional  theory  is  used  to  obtain  the  energy  of  the 
valence  electrons  is  a  spherically  averaged  external  ionic  potential,  and  the 
ionic  repulsive  energy  is  ealeuleted  by  considering  point  lfko  lone  at  fixed 
sites.  In  the  size  range  studied  here  (up  to  90  atoms  per  cluster)  a  layering 
effect  is  observed  in  the  cluster  geometries.  Clutters  are  composed  by  well 
separated  homoatomic  Ne  and  Cs  shells,  the  external  one  being  always  a  Cs 
shell.  There  Is  a  shong  reconstruction  of  the  cluster  geometry  with  increasing 
size.  We  have  also  found  that  the  cohesive  energy  goes  rapidly  to  the  bulk  limit 
and  that  the  mixed  cluster  is  stable  with  respect  to  pure  unmixed  fragments. 
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Electronic  Structure  of  Adsorbates  Near  Surface  Irregularities:  Cluster 
Models  for  CO  on  Cu  Surfaces. 

K.  Hermann,  Abt.  Theorie,  Fritz-Haber-Institut  Berlin,  Germany 

Workfunction  and  FAX  experiments  have  shovn  that  surface  irregularities 
such  as  kinks,  steps,  or  imperfections  can  give  rise  to  large  local  inhomo¬ 
geneous  electric  fields.  Consequences  for  the  chemisorptive  interaction  are 
studied  in  cluster  models  for  CO  adsorbed  on  Cu  surfaces.  The  CO-Cu  inter¬ 
action  is  modelled  by  Cus(l,4)C0  clusters  vith  the  external  electric  field 
simulated  by  that  of  a  point  charge  of  size  q.  Positive  q  values  simulate 
the  situation  near  a  surface  step  where  an  effective  positive  charge  exists 
(Smoluchovski  effect)  vhile  negative  q  simulate  the  field  of  a  negatively 
charged  coadsorbate  ion.  The  electronic  cluster  wavef unctions  are  obtained 
from  ab  initio  Hartree-Pock-LCAO  calculations.  For  positive  q  the  1-elec¬ 
tron  energies  of  all  cluster  orbitals  shift  downwards  which  can  be  ex¬ 
plained  by  a  (dominantly  linear)  surface  Stark  effect.  Further,  the  two  e- 
type  cluster  orbitals  responsible  for  n  backdonation  in  Cu^CO  increase 
their  CO  2rt  admixture  in  the  presence  of  the  field.  The  increased  effective 
CO  2rt  occupation  leads  to  a  field-induced  C-0  bond  weakening.  This  explains 
the  increased  C-0  equilibrium  distance  and  the  reduced  frequency  of  the  C-0 
stretching  vibration  found  in  the  present  model  calculations  which  i9  con¬ 
sistent  with  experimental  findings  for  CO  adsorbed  near  surface  steps.  The 
results  for  negative  q  are  complementary  to  those  for  positive  q  which  is 
consistent  with  the  experimental  result  that  negative  coadsorbate  ions  can 
reduce  the  CO  dissociation  probability  at  metal  surfaces. 
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Hydrogen-Magnesium  Interaction  Energies 

U.  Yxklinten,  Z.  J.  Tian,  A.  Mallo,  C.  Engdahl,  and  B.  I.  Lundqvist 
Institute  of  Theoretical  Physics, 

Chalmers  University  of  Technology ,  S  -  412  96  Goteborg,  Sweden 

Abstract 

Total-energy  calculations  with  the  effective-medium  theory,  EMT,  have 
been  performed  on  the  hydrogen-magnesium  system  in  order  to  obtain 
hydrogen  chemisorption  and  dissolution  properties.  The  EMT  has  been 
applied  to  a  single  H  adatom  on  the  Mg(OOOl)  and  (1 120)  surfaces,  to  a 
layer  of  H  atoms  adsorbed  on  the  Mg(0001)  surface,  to  H  coadsorbed  with 
Mg,  to  an  additional  H  atom  coadsorbed  with  a  full  H  overlaver  on  the 
Mg(0001)  surface,  and  to  H2  on  the  Mg(0001)  surface.  Heat  of  solution  for 
H  in  bulk  Mg.  and  Mg  lattice  expansion  as  a  function  of  the  H  concentration, 
have  also  been  calculated.  At  several  instances  comparisons  with  first- 
principles  calculations  and  experimental  data  can  be  made,  and  the  fair 
agreement  gives  confidence  in  the  method  used.  Out  of  the  calculations 
emerge  a  picture  with  activated  dissociative  chemisorption  of  H;  on  the 
Mg(0001)  surface,  initial  adsorption  in  subsurface  sites,  low  activation 
barriers  for  H-atom  diffusion  on  the  surface  and  in  the  bulk,  and  with  H-H- 
interactions  having  several  important  consequences.  The  calculations  also 
provide  numbers  for  direct  comparison  with  experiment,  such  as. 
adsorbate  vibrational  frequencies  and  diffusion  barriers  for  adsorbed  and 
dissolved  hydrogen,  and  illustrate  the  great  importance  of  allowing  fo: 
relaxation  of  the  atomic  configurations. 
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Ab-initio  calculations  of  relaxations  and  workfunctions  of  transition  met- 
alsurfaces:  The  anomalous  relaxation  of  Pd(100)  and  Rh(100j 

D.  Hennig1,  M.  Methfessel,  S.  Weber,  and  M.  Scheme: 

Fritz-Haber  Institut  der  MPG 
Faradavweg  4-6,  D-1000  Berlin  33,  Germany 

Experimentally,  most  clean  surfaces  show  an  inward  relaxation,  i.e.  the  inte:-iaye: 
spacing  between  the  first  and  second  layer  is  smaller  than  the  bulk  inter-iayer  spacing. 
However,  some  transition-metal  are  found  to  relax  outward.  Current  theories  focus  or. 
the  effect  of  charge-density  smoothing  at  the  surface  and  can  only  explain  inward  re¬ 
laxation.  We  therefore  performed  density-functional  theory  calculations  for  several  fee 
transition-metal  surfaces,  in  order  to  analyse  this  hitherto  non  explained  effect.  In  par¬ 
ticular  we  investigated  the  trend  for  different  low  index  fee  crystal  surfaces:  (111).  (1001. 
and  (110).  We  used  the  full- potential  linear-muffin-tin  orbital  method  and  a  slab  geome¬ 
try.  Calcutated  relaxations  and  workfunctions  agree  well  with  experimental  results  For 
Pd(iOO)  and  Rh(100)  we  find  indeed  the  so-called  anomalous  behavior.  The  results  are 
explained  in  terms  of  a  competition  between  charge-density  smoothing  and  the  bonding 
and  antibonding  nature  of  d-derived  states. 


‘Permanent  address:  Department  of  Theoretical  Physics,  Humboldt-L'niversitv  of  Berlin  In^aiiaenstr 
110.  Berlin  0-1140  Germany 
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Interaction  of  the  Pt(lll)  surface  with  adsorbed  Xe  atoms 
J.E.  Muller 

Institut  fur  Schicht-  urtd  Inrteiuechnik,  Forschungszentrum  Jiilich. 

D-5170  Jiilich.  Federal  Republic  of  Germany 

The  interaction  potential  between  the  Pt(lll)  surface  and  Xe  atoms 
adsorbed  upon  it,  was  calculated  using  the  cluster  method  with  the 
Kohn-Sham  scheme.  For  a  single  Xe  atom,  an  energetically  most  favor¬ 
able  on- top  adsorption  site,  3.0  A  above  a  Pt  atom  was  found,  with  an 
activation  energy  for  surface  diffusion  of  22  meV,  ana  a  perpendicular 
vibration  energy  of  8.5  meV.  This  vibration  energy  is  reduced  drastically 
by  clustering  of  the  adsorbed  Xe  atoms. 
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ADSORPTION  AND  MAGNETISM . ICGTO-MP-LSO  STUDY  OF  ADATOMS  ON  FERROMAGNETIC 

Ni(lOC)  SURFACE. 

Franca  Mele,  N. Russo  and  Marirosa  Toscano 

Dipartimer.to  di  Chimica.Universita  della  Calabria,  I-8703C  Arcavacata  di 

Rende(CS),  Italy. 

The  relationships  between  adsorption  of  adatoms  (H,  C,  N,  0  and  CCi  anc 

changes  in  the  saturation  magnetization  Ms  of  ferromagnetic  Ni ( 100 (surface 

have  been  investigated  with  the  first-principle  LCGTO-MP-LSD  method. 

Our  results  show  that: 

i.  the  most  stable  adsorption  sites  for  H,  C,  N  and  0  atoms  is  the  four¬ 
fold  hollow  one,  while  CO  adsorbs  preferentially  in  on  top  position. 

ii.  the  computed  geometric  and  spectroscopic  equilibrium  parameters  are  in 
good  agreement  with  available  experimental  data. 

iii.  the  spin  polarized  calculations  show  that  H,  C  and  N  adsorption  causes 
a  decrease  of  1 ,  2  and  3  g  in  the  total  magnetic  moment.  In  the  case 
of  0  and  CO  adsorption  no  change  in  the  total  magnetic  moment  is  ob¬ 
served  but  significant  variations  are  found  for  the  local  atomic  mo¬ 
ments  of  surface  Ni  atoms  directly  involved  in  the  interaction  with 
the  adatoms . 
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A  Study  of  Alkali  Adlayars  on  Al{111)  as  a  Function  of  Coverage 

P.  A.  Serena 

IBM  Research  Division,  Zurich  Research  Laboratory 
8803  Riischlikon,  Switzerland 
and  N.  Garcia 

Departamento  de  Fisica  ds  la  Materia  Condensada 
Universidad  Autonoma,  28049  Madrid,  Spain 

We  have  studied  the  adsorption  of  Na,  K  and  Cs  on  Ai(ttt)  as  a  function  of 
coverage  0  by  performing  self-consistent  calculations  on  thick  slabs  where  the 
local  Ionic  pseudopotentials  are  first  averaged  parallel  to  the  surface.  As  0  is 
Increased,  iiie  difference  charge  density  induced  by  the  adlayer  dramatically 
shows  electron  transfer  towards  the  substrate  up  to  a  critical  ©e~0.3  ML, 
followed  by  a  slower  accumulation  in  the  outer  adlayer  region.  The  work  function 
change  A<t>  exhibits  a  concomitant  asymmetric  minimum  (In  agreement  wit’ 
experiments),  So  dona  Alan  fhi*.  muiKhriiim  arllnyer.substrata  distance,  whorooc 
the  frequency  co  of  relative  perpendicular  vibrations  shows  a  maximum  and  is 
quite  low,  especially  for  Cs  This  imp':es  significant  shifts  and  smearing  ot  o  and 
A<t>  at  room  temperature.  Agreement  with  three-dimensional  calculations  for  Na 
adlayers  on  jelllum  can  be  obtained  for  a  reasonable  choice  of  the  alkali  core 
radius  rc.  The  binding  energy  monotonically  decreases  as  ©  is  increased.  It  is 
rather  insensitive  to  rc  and  in  reasonable  agreement  with  the  above-mentioned 
calculations. 
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Ab-initio  calculations  of  chemisorption  and  diffusion  of  AI  and  Si  on  Al(lll) 
R.  Stumpf,  J.  Neugebauer  and  M.  Scheffler, 

Frits-Haber-Institut  der  MPG,  Berlin,  Germany 

We  studied  the  clean  Al (111)  surface  and  the  chemisorption  of  Al  and  Si  atoms  at  low 
and  at  high  coverages  (©  =  0.1-1  ML)  on  this  surface.  We  therefore  used  density- 
functional  theory,  a  plane  wave  basis  set  for  the  electronic  wave  functions,  the  local  density 
approximation  for  the  exchange-correlation  functional  and  ab-initio,  norm-conserving, 
fully  separable  pseudopotentials.  Calculating  self-consistently  the  electronic  structure, 
we  simultaneously  relaxed  the  adsorbate  and  the  two  top  surface  layers.  We  discuss  the 
nature  of  the  bonding,  which  is  mainly  covalent,  and  the  rather  large  atomic  relaxation 
around  the  adsorbate.  The  relaxations  in  our  calculations  are  clearly  more  important  than 
a  recent  article11  about  adsorption  on  the  Al(100)  surface  indicates.  These  unexpected 
results  have  important  consequences  for  the  binding  energy  at  different  sites  and  for 
surface  diffusion. 

*P.  J.Feibelman,  Phys.  Rev.  Lett.  65,  729  (1990) 
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calculation  of  migration  energies  for  DIFFUSION 
ON  THE  Si(IOO)  SURFACE. 

G.  Brocks.  P.J.  Kelly. 

Philips  Research  Laboratories.  5600  JA  Eindhoven.  The  Netherlands. 

R.  Car 

S1SSA,  Trieste,  Italy. 


The  Car-Parrinello  approach  to  electronic  structure  calculations  enables  us  to  use 
supercclls  containing  a  sufficiently  large  number  of  atoms  that  we  can  study  the 
interaction  between  adatoms  and  a  surface  in  the  limit  where  the  adatom-adatom 
interaction  is  negligible.  Results  are  presented  for  the  binding  position  and  the  mi¬ 
gration  harriers  of  single  Si.  P  or  A1  adatoms  on  the  reconstructed  Sit  1 00 1  surface 
These  results  arc  extracted  from  an  energy  surface  which  is  constructed  hv  calculating 
the  total  energy  for  a  number  of  position  -  of  the  adatom  on  the  surface.  At  each  po¬ 
sition.  the  total  energy  is  minimized  simultaneously  with  respect  to  the  electronic 
charge  density  and  all  remaining  (adatom  and  substrate)  ionic  degrees  of  freedom.  Trc 
calculations  are  carried  out  within  the  local  density  approximation,  using  plane  wave* 
and  pseudopotentials  in  a  supercell  conti ining  96  substrate  atoms  plus  the  adatom. 
We  predict  that,  at  room  temperature,  t  te  diffusion  of  P  and  Si  adatoms  is  quasi- 
onc-dimcnsional  in  the  direction  parallel  to  the  dimer  rows  making  up  the  Sit  I  (XT 
surface.  Dimers  of  adatoms  are  found  to  form  stable  structures  on  the  Sit  1 00  j  surface. 
The  binding  energy  and  geometry  of  a  number  of  Si  ad-dimer  structures  are  ealculatcc 
and  their  structural  stability  is  investigated  by  means  of  ab  initio  molecular  dynamic* 
simulations. 
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Atomic  Arrangement  of  Alkali  Adatoms  on  Si(001)-2xl  Surface 

K.  KOBAYASHI,  b.  BLUGEL*,  K.  TERAKURA 

ISSP,  University  cf  Tokyo,  Tckyo  105,  Japan 

♦IFF,  rcrspnungszentruir,  Jiilicn,  D-517C  JOlicr, ,  Germany 

Adsorption  sites  of  Na  aud  Li  on  the  Si(0Gl )-'2x  i  surface  at 
the  saturation  coverage  are  determined  theoretically  by  using  the 
simulated  annealing  method  proposed  by  Williams  and  Sole:.  For 
N‘a.  the  most  stable  configuration  is  the  combination  of  the  hollow 
site  along  the  dimer  chain  and  the  trough  site  on  top  of  the  third  layer 
Si  atom  in  contradiction  with  an  earlier  calculation.  The  substrate 
relaxation  is  a  crucial  factor  for  stabilizing  this  configuration.  Ou 
the  other  hand.  Li  atoms  are  adsorbed  at  low  symmetry  positions. 
For  both  of  Xa  and  Li.  we  observe  a  significant  movement  of  tire 
substrate  atoms  towards  the  idea!  surface  configuration. 


AB  INITIO  MOLECULAR  DYNAMICS  WITH  THE  CONJUGATE  GRADIENT 
METHOD  :  APPLICATION  TO  GERMANIUM 


V.Yu.Milman,  Xudong  Weng,  M.C.Payne,  V.Heine,  R.J.Needs 
Cavendish  Laboratory,  University  of  Cambridge,  Cambridge  CB3  OHE,  UK 

The  use  of  ab  initio  Car-Pamnello  type  molecular  dynamics  (MD)  simulations  is  playing 
an  increasingly  important  role  in  the  understanding  of  dynamical  properties  of  solids  and 
liquids.  In  this  contribution  we  describe  an  implementation  of  the  MD  technique  based  on  the 
conjugate- gradient  method  for  electronic  minimization  [1].  This  method  differs  from  the  original 
Car-Paninello  idea  in  that  it  is  not  necessary  to  introduce  an  arriucial  dynamics  for  the  electronic 
degrees  of  freedom.  Due  to  the  high  numerical  stability  and  rapid  convergence  of  the  conjugate- 
gradient  algorithm  it  was  possible  to  introduce  two  different  time  scales  for  the  MD  simulation, 
performing  the  electronic  minimization  with  a  time  step  which  was  longer  than  the  step  used  for 
the  updating  of  Ewald  forces  on  the  ions. 

Test  MD  runs  were  performed  for  germanium  at  temperatures  close  to  the  melting  point. 
We  used  the  Starkloff-Joannopoulos  local  pseudopotential  and  a  r-point  sampling  for  several 
different  supercells  and  cutoff  energies.  A  Nos£  thermostat  in  Hoover's  formulation  was  used 
to  hold  the  temperature  constant,  and  the  equations  of  motion  of  the  ions  were  integrated  using 
the  Gear  predictor-corrector  method.  Preliminary  results  for  structural  and  thermodynamic 
properties  are  in  reasonable  agreement  with  the  experimental  data  for  cither  liquid  or  amorpnous 
germanium  depending  on  the  temperature  and  the  prehistory  of  the  run. 

We  envisage  subsequent  implementadon  of  this  technique  with  the  transition  state  theory 
to  obtain  free  energies  and  diffusion  coefficients  at  low  temperatures  where  the  characteristic 
time  for  diffusion  is  much  longer  than  the  possible  length  of  an  ah  initio  run. 

1.  M.P.Teter,  M.C.Payne,  D.C.AUan,  Phys.  Rev.  B40  12255  (1989). 
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Application  of  Car-Parrinello  Method  tr^  Solid  and  I. : q •  •  fi 
Metal  and  Metallic  Allov  Guo-Xin  QIAN,  M.  WEINERT,  G.  W.  FER¬ 
NANDO  and  J.  W.  DAVENPORT  Brookhaven  National  Laboratory*1  — 
We  have  performed  first-principles  molecular  dynamics  simulations  for 
solid  and  liquid  metals  such  as  Na  and  Al,  and  for  solid  metallic  alloy 
such  as  Ai-Sn.  We  suggest  a  variable  occupation  scheme  to  soive  th  ■ 
problem  of  degenerate  states  near  fern:  level.  Simulations  were  carrier 
out  at  several  different  temperatures  for  each  metal.  The  calculated 
static  ar.d  dynamic  properties  such  as  pair-distribution  function  g'r- 
coefficient  of  self-diffusion  D  and  activation  energy  Ea  of  liquid  Na  anct 
Al  arc  either  in  good  agreement  with  X-ray  ar.d  neutron  scattering  mea¬ 
surements  or  suggest  further  experiment'  verification.  Furthermore,  we 
studied  a  metallic  alloy  Al-Sti  using  the  same  technique.  We  have  cal 
culated  its  ece:g\  versus  temperature  in  the  solid  phase  which  yields  . 
specific  heat  W  3.5  ;:n  slightly  above  'he  Dulor.tr- Petit  value  of  3kp. 


“Work  suprartia  by  thr  Division  a*  Material'  Zr.tv.iti.  UFDOLi.  ur-.irr 
contract  DE-  AC02-76CH00015. 


Guo-Xin  Qian 
Department  of  Physics 
Brookhaven  National  Laboratory 
Upton,  NY  11973 
U.  ?  A. 
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Electronic  and  Structural  Properties  of  Liquid  Gallium 

X.  G.  Gonga,i,  Guido  L.  Chiarotti®'*,  M.  Parrinello®'*,  E.  Tosatti®4 

a)  International  School  for  Advanced  Studies  (SISSA),  Strada  Costiera  11,  1-34014 
Trieste,  Italy 

b)  International  Centre  for  Theoretical  Physics  (ICTP),  P.0.  Box  586,  1-34014 
Trieste,  Italy 

c)  Laooratorio  Tecnologie  Avanzate  Superfici  e  Catalisi(TASC),  del  Consorzio  In- 
teruniversttar--'  Nazionale  per  la  Fisica  della  Materia  (INFN),  Padriciano  99, 
1-34011  Trieste,  Italy 

i)  IBM  Research  Division,  Zurich  Forsehunglaboratorium,  CH-8803  Ruschlikon, 
Switzerland 

Abstract 

We  have  studied  the  high  temperature  (T  =  ’000k)  properties  of  Liquid  Gallium 
(1-G&),  by  using  the  first  principle  molecular  dynamics  simulation  [1;  at  three  dif¬ 
ferent  densities  (change  of  density  is  chosen  instead  of  temperature,  because  of  the 
shorter  equilibration  time).  The  structure  factor  S(k)  and  pair  correlation  func¬ 
tion  computed  at  the  experimental  density  agree  well  with  available  experiment 
at  T =959K  [2]  •  Study  of  direct  space  configuration  reveals  that  about  half  of 
Ga  atoms  are  instantaneously  engaged  in  a  covalent  bond.  Of  these  again  only 
half  form  binary  bonds,  while  the  rest  form  3-aioms  chains  and  occasionally  4- 
atom  chains  or  4-atom  rpJ  "star”.  This  offer  a  very  instructive  picture  of  how 
1-Ga  mediate  between  the  physics  of  solid  a-Ga  (binary  bond)  and  that  of  d-Ga 
(infinite  chains)  [3].  As  the  density  increase,  a  shoulder  devoleps  in  S(k),  similar 
to  that  observed  experimentally  at  low  temperatures.  We  find  a  direct  connection 
between  the  presence  of  this  shoulder  and  the  proportion  of  chains  over  binary 
bonds  which  we  observed  at  higher  densities.  We  also  provided  a  first  ab-initio 
description  of  other  electronic  properties  of  liquid  Ga. 

Reference 

1.  Car  R.  and  Parrinello  M.,  Phvs  Rev  Let, 55  2471(1985). 

2.  Bellissent-Funel  M.  C-,  Chienx  P.,  Levesue  D.  and  Weis  J.  J.,  Phys  Rev  A  39 
6310(1989). 

3.  Gong  X.G.,  Guido  L.  Chiarotti,  Parrinello  M.,  Tosatti  E.,  To  be  published. 
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LMTO  based  methods  for  local  electronic  structure 
in  substitutional  and  structural  disorder 

T  Jarlborg  (1),  E^Moroni  (2),  J  Haglund  (3)  and  J  Sjostrom  (3) . 

(1)  DPMC,  University  of  Geneva,  Switisrldod 

(2)  IRAMA  and  University  of  Lausanne,  Swiierland 

(3)  Royal  Inst  of  Technology,  Stockholm,  Sweden 


Abstract 

Methods  for  calculations  c l  the  local  electronic  structure  at  a 
given  site  surrounded  by  disorder  of  various  kinds  are  developed. 
The  methods  are  based  on  the  self-consistent  U4T0  method.  In  one 
application  substitutional  disorder  as  in  alloys  is  considered  and 
the  local  electronic  structure  at  each  constituent  is  obtained.  The 
results  are  compared  with  results  frcm  U4T0  super-cell  calculations 
and  with  published  CPA-calculations.  In  a  second  application,  we 
consider  structural  disorder  as  in  amorphous  systems,  in  liquids  or 
in  solids  at  elevated  temperatures.  Results  for  Pd  with  disorder 
around  the  basic  fee  structure  are  presented  and  compared  to  other 
calculations  and  to  experiment.  The  results  include  examples  of 
density-of-states  functions  and  total  energies,  and  the  poasibe  use 
from  further  development  of  the  two  methods,  or  from  a  combination 
of  them,  is  discussed. 
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MICROSCOPIC  APPROACH  TO  THE  ATOMIC  AND  ELECTRONIC  STRUCTURE  OF  TRANSITION-METAL  GLASSES 

•CH.  HAUSLEITNER  AND  W.  JANK 

Institut  fOr  Theoretische  Physik,  TU  Wien,  Wiedne-  HauptstraDe  8-lo,  A-loAo  Wien 

He  present  a  novel  hybridized  NFE-TB8  approach  for  the  interatomic  forces  in 
disordered  transi tion-metal  alloys'to  construct  realistic  models  for  the  ato¬ 
mic  structure  of  amorphous  alloys.  These  pair  potentials  show  clearly  the 
Interrelation  between  short-range-order  (SRO)  and  the  shape  and  filling  of  the 
d-band.  Here  we  present  calculations  of  the  atomic  and  electronic  structures  of 
amorphous  transition-metal -alloys  at  the  example  of  NixZr^_x  (x=o.35,  o.5o,  o.65) 
based  on  molecular  dynamics  simulations,  supercel 1-electronlc  structure  calcu¬ 
lations,  and  on  a  single-scatterer  final-state  approximation  for  the  photo- 
omission  spectra. 
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Effect  of  disordered  Be  addmixture  upon  pre-Fermi  electron 
spectra  and  constant  energy  surfaces  in  high-T  oxide 
La,  Ba  CuO  c 

f -X  x  4 

A.V .Zhalko-T itarenko ,  V . V . He ao shkal enk o ,  S . N. Cher nukh in 
L.H . She  1 udc henk o 

As  it  was  previously  exper i mental  1 y  established  by  Simo¬ 
nov  at  al .  ,  high-Tc  oxide  of  interest  demonstrate  superconduct¬ 
ivity  at  x‘0,12  only  when  admixture  is  disordered  .  Thus  one 
can  not  expect  to  adequately  describe  its  electronic  structure 
via  traditional  band-structure  approach.  We  have  carried  cut 
cluster  calculation  -for  perovskite  sub  cells  of  its  structure, 
containing  both  La  and  Ba  atoes.  Cluster  level  half  -  width 
was  adjusted  to  band-structure  data  for  x*0,  and  was  treated 
as  proportional  to  x  for  xWO.  For  electronic  structure  was 

treated  as  weighted  admixture  of  correspondent  broaden  La-  and 
Ba-  containing  cluster  energy  levels.  Calculation  results 
allowed  us  to  state,  that  for  optimal  x*0,l7  La  Ba  CuO 

electron  spectra  demonstrate  a  sharp  oxygen  peak  \nK a  *clos$ 
vicinity  of  Fermi  level  with  energy  less,  than  phonon  limiting 
frequency.  This  peak  is  dominantly  formed  by  La  layer  oxy¬ 
gen.  Pre-Fermi  constant  energy  surfaces  where  evaluated  via 
Fourier  transf ormation  of  real-space  charge  density.  There 
where  not  found  any  pronounced  changes  in  its  structure  due  to 
disordered  Ba  addition  exactly  at  . 
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